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Along this review, we focus on the study of several properties of modified gravity theories, in 
particular on black-hole solutions and its comparison with those solutions in General Relativity, 
and on Friedmann-Lemaitre-Robertson- Walker metrics. The thermodynamical properties of fourth 
order gravity theories are also a subject of this investigation with special attention on local and global 
stability of paradigmatic f{K) models. In addition, we revise some attempts to extend the Cardy- 
Verlinde formula, including modified gravity, where a relation between entropy bounds is obtained. 
Moreover, a deep study on cosmological singularities, which appear as a real possibility for some 
kind of modified gravity theories, is performed, and the validity of the entropy bounds is studied. 

PACS numbers: 98.80.-k; 04.50.+h 

I. INTRODUCTION 

General Relativity (GR) has been the most successful gravitational theory of the last century, fully accepted as 
a theory that describes the macroscopic geometrical properties of space-time. For an isotropic and homogeneous 
geometry, GR leads to Friedmann equations vifhich describe in an appropriate way the cosmological evolution with 
radiation and then matter dominated epochs. Nevertheless, the development of observational cosmology in the last 
decades with experiments of increasing precision like supernovae observations |IH3] has revealed that the Universe is 
in a stage of accelerated expansion. GR provided with usual matter sources is not able to explain this phenomenon. 
Moreover, GR does not account either for the cosmological era known as inflation believed to have taken place 
before the radiation stage and that could alleviate some problems of standard cosmology like the horizon and the 
flatness problem [5]. In addition, GR with usual baryonic matter cannot explain the observed matter density deter- 
mined by fitting the standard ACDM model to the WMAP7 data [5], the latest measurements from the BAO (Baryon 
Acoustic Oscillations) in the distribution of galaxies [7] and the Hubble constant (i/o) measurement [8^. Thus, GR 
requires the introduction of an extra component called dark matter that accounts for about 20% of the energy content 
of our Universe. 

A more puzzling problem is associated to the present accelerated expansion of the Universe. There are also a 
large amount of different explanations. One of them, assuming the validity of GR, postulates the existence of an 
extra cosmic fluid, the dark energy (DE), whose state equation p — wdeP (where p and p are the pressure and the 
energy density of the fluid) demand cjde < ~l/3 in order to provide an accelerated cosmic expansion [51111|. The 
cosmological constant is the simplest model of DE, corresponding to an equation of state wde = — 1- However, if we 
assume that the cosmological constant represents the quantum vacuum energy, its value seems to be many orders of 
magnitude bigger than the observed one [T2]. 

Thus, alternative explanations to the cosmological constant have been largely studied in the last 
years. These theories essentially modify GR by considering actions different from the Einstein-Hilbert 
one [I3ti24j . Examples are Lovelock theories, free of ghosts and whose field equations contain second derivatives 
of the metric at most; string theory inspired models, which include a Gauss-Bonnet term in the Lagrangian (see |25| 
and references therein); scalar-tensor theories like Brans-Dicke one, in which gravitational interaction is mediated by 
both a scalar field and GR tensor field; or the so-called f{R) theories (see [26H33] for a recent and extensive review). 
In this investigation we shall restrict ourselves to f{R) theories in the metric formalism (where the connection depends 
on the metric, so the present fields in the gravitational sector of the action come only from the metric tensor) in the 
Jordan frame. In this frame, the gravitational Lagrangian is given by i? -I- f{R), where f{R) is an arbitrary function 
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of the scalar curvature R, and Einstein's equations usually become fourth order on the metric derivatives. The f{R) 
theories were proved (see [261 - 138] among others) to be able to mimic the whole cosmological history, from inflation to 
the actual accelerated expansion era. Diverse applications of these theories on gravitation and cosmology have been 
also widely studied [551 - HT] , as well as multiple ways to observationally and experimentally distinguish them from GR. 
Concerning local tests of gravity and other cosmological constraints, see [32H11]- 

In fact, most of the modified gravity theories present the so-called degeneracy problem: from large scale observations 
(la type supernova, BAO, or the cosmic microwave background) which depend uniquely on the evolution history of 
the Universe, the nature and the origin of DE cannot be determined due to the fact that identical evolutions can be 
explained by a diverse number of theories. However, it has been proved in [TflfiS] and more recently in that when 
scalar cosmological perturbations are studied, f{R) theories, even mimicking the standard cosmological expansion, 
provide a different matter power spectrum from that predicted by the ACDM model [50J. 

The study of alternative gravitational theories to GR requires to confirm or discard their validity by obtaining 
solutions that can describe correctly, e.g., the cosmological evolution, the growth factor of cosmological perturbations 
and the existence of GR-predicted astrophysical objects such as black holes (BH). Therefore, it is interesting to study 
the properties of BH in this kind of theories, since some of their features might be either exclusive of Einstein's 
gravity or intrinsic features of any covariant gravitational theory. On the other hand, obtained results could provide a 
method to discard models that disagree with expected physical results. In this sense, research of BH thermodynamics 
may shed some light about the viability of alternative gravity theories since local and global stability regions, and 
consequently the existence itself of BH, depend on the values of the parameters of the model under consideration. 

BH properties have been widely studied in other modified gravity theories: for instance [FTJ [S2] stud- 
ied BH in Einstein's theory with a Gauss-Bonnet term and a cosmological constant. Gauss-Bonnet 
and/or quadratic Riemann interaction terms are studied in |53) . where it is found that for a neg- 
ative curvature of the horizon phase transitions might occur. BH in Lovelock gravitational theories 
were studied in [541 155) . where the corresponding entropy was calculated. Other recent works have 
studied [56l - f58] BH in the context of Hofava-Lifshitz gravity and static solutions in the context of /(T) theories 
[55] as well. 

Previous works concerning BH in f{R) theories proved [5D that for a Lagrangian R + aR^ the only spherically 
symmetric solution is Schwarzschild's one provided that one works in the Einstein's frame. Again in Einstein's frame, 
[61] proposed uniqueness theorems for spherically symmetric solutions with an arbitrary number of dimensions (see 
[52] for additional results). 

Spherical solution with sources were also studied in [63] whereas |64j developed a new covariant 
formalism to treat spherically symmetric space-times claiming that Schwarzschild solution is not a unique static 
spherically symmetric solution. Spherically symmetric /(i?)-Maxwell and /(_R)-Yang-Mills BH were studied in [65], 
confirming the existence of numerical asymptotic solution for the second ones. Concerning axially symmetric solutions, 
authors in |66j showed that these solutions can be derived by generalizing Newman and Janis method to f{R) theories. 
An scalar-tensor approach is used in [67] to show that Kerr BH are unstable in a subset of f{R) models because of 
the superradiant instability. In [68" the entropy of BH is calculated in the Palatini formalism by using the Noether 
charge approach. Anti-de Sitter (AdS) BH have been studied [69j in f{R) models using the Euclidean action method 
(see, e.g., [7D1[7T]) to determine different thermodynamic quantities. In [72], the entropy of Schwarzschild-de Sitter 
(SdS) BH is calculated in vacuum for certain cosniologically viable models, and their stability is discussed. 

Different aspects of BH in /(i?) have been studied recently, including in particular stability issues [65, 73 , showing 
an increasing interest in these topics, static and spherically symmetric solutions [73H75] . rotating configurations [771175] 
and the existence of anomalies 179] . 

In fact, since the establishment of thcrmodynamical properties for black holes along the seventies, a special interest 
has been focused on the relation of the Einstein gravity and laws of thermodynamics, whose efforts have been increasing 
along recent years. Particularly, Jacobson showed in 1995 ^80j that field equations in General Relativity can be 
recovered from thermodynamics identities by assuming the relation, given by mechanical laws of black holes, between 
the entropy and the area of the horizon, and assuming that the relation 5Q = TdS holds for all local Rindler 
causal horizons through each spacetime point. Such an important result points to the hypothesis that gravitational 
equations may be just equations of state of the spacetime rather than fundamental laws of gravitation. Furthermore, 
this suggestion has been reinforced since the relationship between gravitational equations and thcrmodynamical ones 
has been extended to more general gravitational theories, as for instance, f{R) gravity (see [STj). 

Moreover, Friedmann-Lemaitre-Robertson- Walker (FLRW) equations have been successfully derived from the ther- 
modynamical properties of the apparent horizon (see [52HM] ). In this sense, modified FLRW equations of more general 
gravitational theories as Lovelock gravity or F(R) gravity have been also related with the thcrmodynamical proper- 
ties of the apparent horizon !85r'96j. In addition, for a radiation dominated universe, the FLRW equations can be 
rewritten as an entropy relation, which gives an equivalence with the Cardy formula from a two dimensional conformal 
field theory (CFT) (see [57]). This relation, proposed initially by Verlinde in [55], assumes an entropy bound for the 
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universe related with the entropy proportional to the area of the size of the universe, a kind of holographic princi- 
ple. Finally, the formula may be rewritten as a dynamical entropy bound from which a number of entropy bounds, 
proposed earlier, follow. However, such a relation with the 2d CFT has been attempted to be extended for general 
perfect fluids (see |99j), and more general scenarios including modified gravity flOO], but the equivalence could not 
be extended. Furthermore, it was found that the bound on the entropy may be violated in some particular scenarios, 
mainly when the universe moves close to a future singularity. Moreover, in spite of that in general the bound is 
violated sufficiently close to the singularity, even when some quantum effects are included, the violation can be also 
produced much before the singularity occurs [100] . 

Hence, all of these suggest a deep connection between gravitational physics and thermodynamics, which is still an 
open issue nowadays, and could help to improve our knowledge on the way to reconstruct consistent gravitational 
theories that may solve some of the most important problems on theoretical physics. Along this work, some of these 
results are reviewed, and we may suggest new scenarios where gravitational physics may be tested 

The present review is organized as follows: In Section |ll] we revised the rudiments of the f{R) modified gravity 
theories that will be used throughout the paper. Then, Section|lIl]is devoted to study static and spherically symmetric 
configurations in these theories, with special focus on electromagnetic solutions and perturbative approach. The fol- 
lowing Section IV deals with the Kerr-Newman black holes in f{R) theories. Section|v]encompasses a general revision 
of thermodynamics of f{R) black holes for both anti-de Sitter and Kerr-Newman scenarios. The therniodynamical 
stability is then studied for two paradigmatic examples of f{R) models in Section |VI[ Section VII deals with the 
reconstruction of F{R) actions for some particular cosmological solutions and how this kind of theories are capable 
to reproduce the whole cosmological evolution. In Section VIII we review some of the thermodynamical properties 
of FLRW equations, and its derivation from the first law of thermodynamics, which can be easily extended to F{R) 
gravity. Section |IX| is devoted to some extensions of the Cardy-Verlinde formula to more complex configurations of 
FLRW universes, including modified gravity. In Section [Xj future singularities are studied in the context of entropy 
bounds, and the possible violation of the proposed bounds on the energy and entropy of the universe. Finally, we 



conclude the paper by giving our conclusions in Section XI 



II. F{R) THEORIES OF GRAVITY 

In this section, F{R) gravity is introduced, which basically consists in an extension of the Einstein-Hilbert action by 
considering more complex gravitational actions consisting on functions of the Ricci scalar, and which can be rewritten 
in terms of a Brans-Dicke-like theory as shown below. The action for i^(i?)-gravity is given by: 

d^x^{FiR) + 2K,^C^) (!) 

Here Cm denotes the Lagrangian of some kind of matter, and = SttG is the gravitational coupling constant. In the 
present work we employ the natural units system in which h = c = 1. Note also that our definition for the Riemann 
tensor is: 

Field equations are obtained by varying the action Equation ([T]) with respect to 5^1/ , 

Ry..F\R) - \g^uF{R) + g^,UF'{R) ~ V^V,F'(i?) = k'T^J^^ (3) 

where T^i™'' = — y=|^n7 is the energy-momentum tensor and prime denotes here derivatives with respect to R. 
Nevertheless, in general it is very difficult to get exact solutions directly from field Equation (|3]), although some 
particular solutions can be obtained, as shown in the next sections. However, the action Equation ([T]) is equivalent 
to a kind of scalar-tensor theory with a null kinetic term, which can be used to reconstruct the appropriate action for 
some particular solutions, specially for cosmological solutions, as shown in Section VII Then, the action Equation ([T]) 
is rewritten as follows |101| : 



5 = y A^x^g {P{^)R + Q(0) -f 2K^Lm) (4) 
By the variation of the action with respect to the metric tensor g^j^, the field equation is obtained: 

- \9^. {P{4>)R + Qm + P{m>.u + g^..nP{<j>) - V^V,P(0) = K^T^™) (5) 
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In addition, the action Equation Q gives an extra equation for the scalar field 0, obtained directly from the action 
by varying it with respect to 4>: 

P'(0)i? + Q'(0) = (6) 

where the primes denote derivatives with respect to <f>. This equation can be solved such that the scalar field is a 
function of the Ricci scalar R, cj) — 4>{R), and then, by replacing this result in the action Equation Q, the action 
Equation ([T]) is recovered: 

F{R)^P{c^{R))R + Q(cj){R)) (7) 

With the aim of proposing a realistic alternative to GR, a possible modification consists of adding a function of the 
scalar curvature, f{R) to the Einstein-Hilbert (EH) Lagrangian. Therefore the gravitational sector of the action 
Equation ([l]) becomes: 



Sg^ / d^xV|5|(i? + /(i?)) (8) 



(1 + fR)R^.u -7^{R + f{R))g^. + ■D^.h = SttGT^, (9) 



By performing variations with respect to the metric, the modified Einstein Equation (|3| turns out to be 

1 
2 

where fn = df{R)/dR and 2?^^ = g^jyD — V^Vy with □ = Vq,V" and V is the usual covariant derivative. 

These equations may be written a la Einstein by isolating on the l.h.s. the Einstein tensor and the /(i?) contribution 
on the r.h.s. as follows |102| : 



SttGT^, - V^Jr + i (/(i?) - RJr) g^. 



(10) 



We can also find the expression for the scalar curvature by contracting Equation ([9| with g^^ , which gives: 

(/;,-l)i?-2/(i?)+3n/fl =8^Gr (11) 

Note that, unlike GR where R and T are related algebraically, for a general /(i?) those two quantities are dynamically 
related. Thus, in principle, non-null curvature solutions are possible in vacuum for certain /(i?) models. 

In addition, one may study the possible corrections at local scales of the Newtonian law and the violations of the 
equivalence principle by using again an auxiliary scalar field A and applying a conformal transformation. Then, we 
may rewrite the action Equation ([s]) as: 



Sg= / d4xV^[(l + /'(A))(i?-A)+A + /(A)] (12) 



Then, by varying this action with respect to A, it gets R — A, and the action Equation ( jsj) i s recovered. By the 
conformal transformation g^^y = e'^g^^, where a = — ln(l + f'{A)), the action Equation (12) is rewritten in the 
so-called Einstein frame: 



S§ = I d*^ 



R-^d^ad^a~V{a) 



(13) 



where the tilde represents that the scalar curvature is evaluated with respect to the transformed metric 5^^, and the 
potential V{a) is given by: 

M'T) A + fjAja)) 

^ ^ l + (l + /'(A(a)))2 ^ ) 

Here, one has to solve the equation a — — ln(l + f'{A)) in order to get A — A{a). This conformal transformation 
affects the matter sector in the gravitational action by a term proportional to e^°^, which can imply large corrections 
to the Newtonian law unless the potential Equation ( 14 1 is constructed in order to avoid them at local scales by the 
so-called chameleon mechanism |103l I104j , which implies that the mass of the scalar field a should be enough large 
at local scales in order to avoid large corrections to the Newtonian law |105j : 

^ 1 dV(a) ^ 1 r A{a) A + f{A{a)) 1 \ 

- 2 d(72 2\l + r{A{a)) (l + /'(A(a)))2 + /"(A(a))/ ^ 

Hence, if the mass is much larger at local scales than the curvature, the possible corrections to the Newtonian law 
will be negligible. In this sense, several models that accomplish this rule have been proposed (see [106H110] ). which 
are also capable to reproduce the cosmological history. 
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III. STATIC AND SPHERICALLY SYMMETRIC BLACK HOLES IN f{R) GRAVITIES 

Let us study in this section the f{R) sohitions under staticity and spherical symmetry assumptions. The most gen- 
eral D > A dimensional metric external to static and spherically symmetric — therefore non-rotating — configurations 
can be written as (see 

ds^ ^ e-2*WA(r)di2 - A~\r)dr^ - rMnl_^ (16) 

or alternatively: 

ds^ = A(r)d<2 - ^-i(r)dr2 - r^dnl_2 (17) 

where dr2|,_2 is the metric on the S^^^ sphere and identification A(r) = e^^*'-'"-'v4(r) and fj,{r) — A{r) can be 
straightforwardly established if required. Since the metric is static, the scalar curvature R in D dimensions depends 



only on r and it is given, for the metric parametrization Equation (16), by 



R{r) = -^[D^ -5D + 6 + rA'{r){~2D + 3r^'ir)+4) 

- r^A"{r) - A{r) {D^ - 5D + 2r^<^>'{r f - 2{D - 2)r$'(r) - 2r^<^>"{r) + 6)] (18) 
where the prime denotes derivative with respect to r coordinate. 

A. Spherically Symmetric and Static Constant Curvature Solutions: Generalities 

At this stage it is interesting to ask about which are the most general static and spherically symmet- 



ric metrics with constant scalar curvature Rq. This curvature can be found solving the Equation (18) 
R = i?o- Thus, it is straightforward to check that, for a constant <i>(r) — $o, the general solution is [112[ 1113] : 

A(r) = 1 + ^ + ^ fl9) 

with oi and 02 being arbitrary integration constants. In fact, for the particular case D = A, Rq = and $0 = 0, the 
metric can be written exclusively in terms of the function: 

= 1 + 7 + 5 (20) 

By establishing the identifications oi = —2GnM and 02 — Q^, this solution corresponds to a 
Reissner-Nordstrom solution, i.e., a charged massive BH solution with mass M and charge Q. Further comments 
about this result will be made in Section fill Bl 

B. Spherically Symmetric and Static Constant Curvature Solutions 

Let us study here which are the vacuum constant curvature solutions for f{R) theories provided that the metric under 



consideration is of the form Equation ( 16 ) — or alternatively Equation ( 17 ) . By inserting the metric Equation ( 16 1 into 
the general f{R) gravitational action Sg in Equation ([8|, and making variations with respect to the metric functions, 
A{r) and $(r), the equations of motion become: 

(2 - D){1 + .nR)W{r) - r [f"'{R)R'{rf + f" {R){'^' {r)R' {r) + R" [r))] = (21) 

and 

2rA{r)f"'{R)R'[rf + f" {R)[2DA{r)R' [r) - AA{r)R!{r) + 2rA{r)R"{r) + A'{r)rR!{r)] 
-r{R + f{R)) + g'{R)[-2rA{r)<^'{rf + 2DA{r)^'{r) - 4A(r)$'(r) - rA"{r) + 2rA{r)^"{r) 
+A'{r){2- D + ?.r^'{r))] = (22) 



or alternatively in by using the metric parametrization given by Equation (17) the equations of motion now are 
expressed as: 
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A(r)(l + f'iR)) {2M(r) [{D - 2)A'(r) + rA"(r)] + rX' ir)f,' (r)} 
-2A(r)2 [2^^ir)[{D - 2)i?'(r)/"(i?) + rf(^HR)R' (r)' + rR"{r)f"{R)] 
-r/i(r)A'(r)2(l + f [R)) + 2r\{rf{R + }{R)) = 



rR 



'(r)M'(r)/"(i?)} 



(23) 



and 



A(r)//(r) [2(i^ - 2)A(r) + rA'(r)] (1 + /'(i?)) - 2rA(r)2(i? + f[R)) 

^Ji{r)[2\{r)R'{r) [2{D - 2)A(r) + rA'(r)] /"(i?) + r(l + f'{R)){\'{rf - 2A(r)A"(r))} = 



(24) 



where /', /" and /"' denote derivatives of f{R) with respect to the curvature R. In the foUowing calculations in this 
section, we use the set of Equations ( 21 ) and ( 22 1. In order to simplify the difficulty of the above equations, let us first 
consider here the case of constant scalar curvature R — Rq solutions. Then the equations of motion Equations (21) 
and (22 1 reduce to: 



(2-Z?)(l + /'(i?))$'(r) = 



(25) 



and 



R + f(R) + {l + f'{R)) 



'A"ir) + iD-2)^-i2D-4)^^^ 



-3A'(r)$'(r) + 2A(r)$'2(r)-2A(r)$"(r)] = 
As commented in Section |llj the constant curvature solutions of f{R) theories in vacuum satisfy: 

Df{Ro) 



Ro 



2{l + f'{Ro))-D 



(26) 



(27) 



provided that 2(l + /'(i?o)) 7^ D. Thus from Equation (25) one concludes (excluding pathologic cases) that <I>'(r) = 
and then Equation (|26|) becomes: 



A"{r) + {D-2)^ = -^Ro 
r D 



(28) 



This is a /(i?)-independent, linear second order inhomogeneous differential equation. This statement must be under- 
stood in the sense that Equation (28) does not involve explicitly any /(i?) parameter. Nonetheless, it is obvious from 



Equation ( [27| that the value of Rq does depend upon the f(R) model under study. Therefore the metric solution 
A{r) will also inherit information about the model under consideration. Then, Equation (28) can be easily integrated 
to give the general solution: 



A{r) = Ci + 



Co 



(29) 



D{D-l) 

which depends on two arbitrary constants Ci and C2. However, two important remarks need to be done at this stage: 



• Firstly, by comparison with Equation (19), one can see that the term with the power ^ is absent. This fact 
will be studied in Section IIII CI 

• Secondly, this solution has no constant curvature in the general case since, as we found above, the constant 
curvature requirement demands Ci = 1. This issue just requires a constant fixing (or equivalently a time 
reparametrization) and does not affect the solution. 



Then, for negative i?o, the found solution in Equation ( |29[ ) with Ci = 1 is basically the Z?-dimensional generalization 
obtained by Witten [7T] of the BH in AdS space-time solution considered by Hawking and Page [70] . With the natural 
choice $0 = the solution can be written as: 



A(r) = 1 



R 



D-3 



r 



(30) 
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where 



£,_3 WttGdM 2tt 



D-1 



=00^2)^ ' ^^-2-F7i5^ (31) 



2 



where /i£i_2 is the area of the D — 2 sphere, P = —D{D — l)/i?o is the asymptotic AdS space scale squared and M 
is the mass parameter usuahy found in the hterature. 

Thus, the main conclusion of this analysis |112l I113| is that the only static and spherically symmetric vacuum 
solutions with constant (negative) curvature of any f{R) gravity is just the Hawking-Page BH in AdS space. However 
this kind of solution is not the most general static and spherically symmetric metric with constant curvature, as can be 



seen by comparison with the solutions found in Equation (19 1. Therefore there exist constant curvature BH solutions 
that cannot be obtained as vacuum solutions of any f{R) model. In particular, the term with r^~^ power cannot be 



mimicked by any f{R) model. As we shall show immediately, the most general solution as given by Equation (19) 
can be described as a charged BH solution in a D = 4 /(i?)-Maxwell theory. 



C. f{R) Solutions Combined with Electromagnetism 



Let us briefly mention here for c omp leteness with the previous results the case of charged BH in f{R) theories that 

We shall limit ourselves to the D — A case, since otherwise {D ^ 4) the r^^^ 



we shall study in detail in Section IV 



term appearing in Equation ( 19 1 cannot be originated by the Maxwell action to be presented below. Let us consider 



here a /(i?) gravitational Lagrangian supplemented with the usual Maxwell action for standard Electromagnetism. 
Therefore, the considered action represents a generalization of the Einstein-Maxwell action: 



where as usual Fi, 



g\{R + f{R)-F^,Fn 



(32) 



Considering an electromagnetic potential of the form 



= {y{r),Q) and the static spherically symmetric metric Equation (17), one finds that the solution with con 
stant curvature i?o reads [1121 1113| : 

V{r) 



Q 

r 



\{r) = ii{r) = 1 



2GnM 



1 + f'{Ro) r2 



i?0 2 

12 



(33) 



Notice that unlike the EH case, the contribution of the BH charge to the metric tensor is corrected by a (1 + /'(i?o))^^ 
factor. This fact is of particular importance in order to constrain f{R) models as astrophysically viable: the sign 
of the factor 1 + /'(i?o) needs to be positive in order to guarantee that the solution Equation (33) keeps the same 
sign as in the Reissner-Nordstrom solution. Otherwise, for a charge Q the third term in the solution Equation (33) 



will be negative and therefore this fact becomes experimentally observable (affecting for instance to the geodesic 
movements around a charged BH). In fact, the requirement 1 + f'{Ro) > has been usually accepted ^,1141 as one of 
the gravitational viability conditions for f{R) models and is also related, as will be shown here, with thermodynamical 
viability. 



D. Perturbations Around Schwarzschild— (anti)-de Sitter Solutions 

In the previous section we have revised the static spherically symmetric solutions with constant curvature. In EH 



theory, the constant curvature ansatz Equation ( 27 ) provides the most general static and spherical symmetric solution 
straightforwardly. Without entering into details, this is the well-known Birkhoff-Jebsen theorem result |115|, 1116] . 
However, it is not transparent for this to also be the case in /(i?) theories, i.e., that the most general static and 
spherically symmetric solution for /(i?) theories possesses constant curvature. Furthermore, it seems that a generic 
solution would behave differently in the Jordan and Einstein frames, both related by a conformal transformation, at 
least in the perturbative analysis (see |117|. I118j ) . 

In this section we revise the basic rudiments of the perturbative analysis of the problem originally studied in 
[1121 Ill5] . only in the Jordan frame. In that analysis, it was assumed that the modified gravitational Lagrangian is 
a small perturbation around the EH Lagrangian. Therefore, the arbitrary f{R) function is of the form: 



f{R) = -{D-2)AD + ag{R) 



(34) 



where a <C 1 is a dimensionless parameter and g{R) is assumed to be analytic in a. This last assumption is based on 
the fact that the scalar curvature R is also a-dependent as will be seen in expressions Equation (35). Therefore g{R) 



may also be expanded in powers of a. By assumption and for the sake of simplicity, g{R) is assumed analytic in the 
expansion in powers of a. 

Now, assuming that the A(r) and ^{r) functions appearing in the metric Equation ( |17[ ) are also analytical in a, 
they can be written as follows: 



A(r) 



(35) 



where {Ao(r), /J.o(?')} are the unperturbed solutions for the EH action with cosmological constant, i.e., S(A)dS solution, 
given by: 



1 



rD-3 



Al 



(D-l) 



Ao(r) 



-C2{D ~ 2){D - 1) Mr) 



(36) 



which are the standard BH solutions in a ZJ-dimensional Schwarzschild-(A)dS space-time. The factor C2 can be 
chosen by performing a coordinate t reparametrization so that both functions in Equation (36) can be identified. For 
the moment, the background solutions are kept as given in Equation (36) since the possibility of getting A(r) = fi{r) 
in the perturbative expansion will be discussed later on. 

The perturbative procedure is as follows: by inserting Equations (34 1 and (35) into the general Einstein modified 



Equations (23) and (24), equations for each order in a parameter powers are obtained. Please refer to |1121 1113] for 
details. 

Notice that from the obtained results up to second order in a, the corresponding metric has constant scalar 
curvature for any value of the integration constants. As a matter of fact, this metric is nothing but the standard 
Schwarzschild-(A)dS geometry, and can be easily recast in the usual form by making a trivial reparametrization in 
the time coordinate. 

Therefore, this analysis proved that at least up to second order, the only static, spherically symmetric solutions 
which are analytical in a are the standard Schwarzschild-AdS space-times irrespectively of the f{R) model under 
consideration. This result is valid provided that both the solutions and the gravitational Lagrangian are analytic in 
the a parameter and that g{R) corresponds to a small deviation from the usual EH gravitational Lagrangian. Further 
orders in a"^'^' " can be obtained by inserting previous results in the equations of order 3,4, ... and so on. This way 
one may get the solutions {A3.4....(r),/i3.4^...(r)} with the added difficulty that the corresponding equations become 
more and more complicated. 



To conclude Sections |III| and HID we can summarize by saying that in the context of f{R) gravities the only 
spherically symmetric and static solutions of negative constant curvature are the standard BH in AdS space. The 
same result applies in the general case (without imposing constant curvature) in perturbation theory up to second 
order. However, the possibility of having static and spherically symmetric solutions with non-constant curvature 
cannot be excluded in the case of /(i?) functions, which are not analytical in a parameter. 



IV. KERR-NEWMAN BLACK HOLES IN f{R) THEORIES 



In this section, we shall focus our attention in obtaining constant curvature Rq vacuum solutions for fields generated 
by massive charged objects in the frame of f{R) gravity theories. Hence, the appropriate action to derive the field 
equations will be now Equation ( 32 ) and D = 4 will be assumed throughout this section. The obtained equations (in 

Gn — 



h = ks = I imits) are: 



1 



R^.. (1 + f'{Ro)) - 2 (^0 + f{Ro)) 



F 



1 



;gfii^-faf3 







(37) 



In the studied case = 4, the trace of the previous equation leads again to expression Equation ( 27 ) due to the 
conformal character of the Faraday tensor, i.e., F^^ ^ = 0. 

For these f{R) field equations the axisymmetric, stationary and constant curvature Rq solution that describes a 
BH with mass, electric charge and angular momentum is analogous to the one found by Carter and published for the 
first time in 1973 fll9'. In Boyer-Lindquist coordinates, the metric takes the form |120l 1121] : 



9 



with 



Ro 
12 



2Mr 



(l + /'(i?o)) 



:= + cos^ ( 



Afl := 1 



12 



S := 1 



2 

12" 



(39) 



where M, a and Q denote the mass, spin and electric charge parameters respectively. Notice that, unhke in the GR 
case, the contribution of the charge of the BH to the metric is corrected by a (1 + /'(i?o)) factor. This feature was 
aheady obtained for Reissner-Nordstrom BH presented in expression Equation pS] ) and originally in [1121 1113j . On 
the other hand, the required potential vector and electromagnetic field tensor in Equation (37) solutions for metric 
Equation ( 38 ) are respectively: 



A 



F 



Qr fdt 

— I — — a sm" ' 



2Qr a cos 9 sin 9 , „ 
A dr : d9 A 



dt , , n, d( 
a — -{r + a^) — 



(40) 



To lighten notation, from now on let us use the notation = / (1 + /'{Rq)) as a normalized electric charge 
parameter of the BH. 



A. Event Horizons 



Let us now briefly study the horizon structure of these BH: according to the event horizon definition, i.e., g^^ = 0, 
these hypersurfaces are found as the roots of the equation A^ = 0, i.e., 



12 

Rn 



r 



. 24 M 12 , , -o. 



(41) 



This is an algebraic fourth order equation that, following |120l 1121] . can be rewritten as: 

(r - r_)(r - rj„t)(r - rext)ir - rcosm) = 
where 

• r_ is always a negative solution with no physical meaning, 

• Tint and re,j;t are the interior and exterior horizon respectively, and 

• fcosm represents — provided that it arises — the cosmological event horizon for observers between r^xt and 



(42) 



The latter divides the accessible region from the hidden region for an exterior observer. The Equation (42 1 can be 



solved by using L. Ferrari's method for quartic equations. The existence of real solutions for this equation is given by 
a factor h, which we shall denote as horizon parameter. 




12 (a^ + q2 



18 Rr 



(43) 



In Figure [T] we show the zeros of the fourth order polynomial A^ that lead to different astrophysical configurations, 
depending on the values and signs of Rq and h. 

Let us stress at this stage that, from a certain positive value of the curvature Rq"*^ onward, the h factor goes to 
zero for two values of a. They can be understood as follows: 
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Figure 1: Graphics showing horizons positions as solutions of the equation Ar — taken from [1201 1121] . On the left panel 
(-Ro < 0) the presented cases are: /i > (I), BH with well-defined horizons, dashed with dots), h — (II), extremal BH, 
continuous line) and h < (III), naked singularity, dashed). On the right panel {Rq > 0) the represented cases are: /i < (I, 
BH with well-defined horizons, dashed with dots), h = (II, extremal BH and III, extremal marginal BH, continuous line), 
and h > (IV, naked singularity and V, naked marginal singularity, dashed) . 



• Upper spin bound, a = a„iax for which the BH turns extremal — the interior and exterior horizons have merged 
into a single horizon with a null surface gravity. This is the usual configuration for the BH to become extremal. 

• Lower spin bound, a = a„iin, below which the BH turns into a marginal extremal BH. This value can be 
understood as the cosmological limit for which a BH preserves its exterior horizon without being "torn apart" 
due to the relative recession speed between two radially separated points induced by the cosmic expansion 

In Figure [2] we revisit the original plot shown in |1201I121] . where for certain values of the electric charge parameter 
Q the allowed ranges of the spin a parameter values are presented. 





-0.5 0.0 



Figure 2: The shaded regions, delimited by the upper amax and lower a-min curves, represent the values of a/M for which the 
existence of BH is possible once Ro value is fixed. Panels show for Q/ M = (left) and Q / M = 0.75 (right) on the left 
and right panels respectively. Note that _Ro has dimensions of [length] when normalizing. Original plots at [1201 1121] . 
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V. BLACK HOLES THERMODYNAMICS IN f{R) THEORIES 



As mentioned in the Introduction, the study of BH thermodynamics has drawn a lot of attention when studying 
these configurations in alternative theories of gravity. We shall divide this section in two main parts: the study 
of thermodynamics in AdS case and then the study of KN case as the natural generalization of the former. This 
approach will enable to introduce first the rudiments of the Euclidean method so the reader can familiarize with it, 
and then extend the formulation for more general objects such as the KN BH configurations. Finally, we will provide 
some examples by considering f{R) models that illustrate the procedure for both configurations. In the AdS case, the 
dimension of the space-time will be arbitrary and then fixed to several values in order to illustrate the found results. 
On the contrary, for the KN study, dimension will be fixed from the very beginning to D = 4 for the sake of simplicity. 



A. BH Thermodynamics for AdS Configuration 



In order to study the different thermodynamic quantities for the f{R) BH in AdS space, the temperature is usually 
considered as the basic quantity from which the rest of the thermodynamic quantities are derived. In principle, there 
are two different approaches of introducing this quantity for the AdS solutions. Firstly, one of the possible definitions 
is the one provided by the Euclidean quantum gravity |122H125] that was revised in [1121 1113) yielding: 



Te = -^e-*(''-)A'(rH) (44) 
47r 

Another possible definition of temperature was firstly proposed in jl26| stating that temperature can be given 
in terms of the horizon gravity K. as Tjc = j^. Then, by reminding the metric definition Equation ( [16^ , it is 
straightforward to obtain that Tjc — Te- In conclusion, both definitions lead to the same result for this kind of 
solution. Notice also that in any case the temperature depends only on the behavior of the metric near the horizon 
and it is independent of the underlying gravitational action. This last statement can be understood by remarking 



that either different actions or different field equations giving rise to the same solutions of the form Equation ( 16 ) 
will also present the same temperature. This is not the case for other thermodynamic quantities as will be seen later. 
Taking into account the results in previous sections, the usual analyses for /(i?) theories [11211113] have focused on the 



constant curvature AdS BH solutions with $ = as a natural choice, and A{r) as given by expression Equation (30 1 
Then, both definitions of temperature lead to 



1 _ AttPth 
T ^ (D - l)r|j + (£> - 3)/2 



P^^ = Jl^^^jr^^h^^^ (45) 



Notice that for a given dimension D, the temperature is a function of th only, i.e., it depends only on the BH size. 
In the limit r^j going to zero the temperature diverges as T ~ 1/ th and for ru going to infinite T grows linearly with 
Th- Consequently T has a minimum corresponding to a temperature: 



JiD- l)(L>-3) , , 

To ^ ^^l^^^ 1 (46) 

The existence of this minimum was established in [70] for D — A hy Hawking and Page. More recently, Wittcn 
extended this result to higher dimensions |71j . This minimum will be relevant in order to set the regions with 
different thermodynamic behaviors and stability properties. 



Explicit solutions for A(r^f) = in Equation (30 1 can be found for D = 4,5 and were found 
in [miTT^ . 

Another important temperature value is the one corresponding to a horizon radius rn — I- Let us denote it as Ti 



and from Equation ( 45 ) it is given by 



Notice that for D > 2 we have Tq < Ti . 

In order to compute the remaining thermodynamic quantities, the standard procedure is to consider the D- 
dimensional Euclidean action defined as follows: 

SE^-j^^Jd^'x^iR + fm (48) 
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When the previous expression is evaluated on some metric with a periodic Euchdean time of period /3, it equals /3 
times the free energy F associated to this metric. The computation by Hawking and Page [70], generalized to higher 
dimensions by Witten [7T], was extended to f{R) theories in [75] and |112U113| . In these references the difference of 
this action, when it is evaluated on the BH and on the AdS metric, was computed leading to: 



ASe 



(flo + /(i?o))/3A^g_2 
32Tr{D-l)GD 



(49) 



where Rq = —D{D — 1)/Z^. From this expression it is straightforward to obtain the free energy as F = ASe/P- 
One can also see that provided that the condition —{Rq + /{Rq)) > is fulfilled — which is the usual case in EH 
gravity — we have > for r^f < Z and F < ioi rn > I- 

Once that both temperature and free energy have been calculated, the total thermodynamical energy may now be 
obtained as: 



E = 



dASp 



dp 



(Rq + f{Ro))Mf 
2{D-1) 



(50) 



where M is the mass defined in Equation (31). This is one of the possible definitions for the BH 
energy for f{R) theories but other definitions are available in the literature, see for instance [127] 
for a more general discussion. For the EH action with non-vanishing cosmological constant, we have 
f{R) = —{D — 2)A£) and then it is immediate to find E = M . However this is not the case for general f{R) gravita- 
tional Lagrangians as seen in Equation (50). Notice that positive energy in AdS space-time requires i?o + f{Ro) < 0. 
Now the entropy S can be obtained from the well-known relation S = /3E — PF. Then one gets: 



S 



{Ro + f{Ro))l^An_2irH) 
8iD-l)GD 



(51) 



where ^_d-2 (''//) is the horizon area given by ^_D-2(^ff) = Notice that once again positive entropy 

requires i?o + /(-Rq) < 0. For the EH action non-vanishing cosmological constant, we have i?o + f{Ro) = — 2(D — 1)/P 
and then the famous Hawking-Bekenstein result [128' is recovered: 



5 = 



AD-2irH) 



4G 



D 



Finally we can compute the heat capacity C which can be written as: 

dE dE dvH 



C = 



dT dr„ dT 



Then it is easy to find: 



C = 



-{Rq + .nRo)){D - 2)^lD^2rT'P (D - l)r|, + (D - 3)/^ 
8Gd{D-1) {D - l)rj, ~ {D - 3)P 

For the particular case of the EH action with cosmological constant we find: 

(i^ - 2)Mp_2rg-' (D - + (D - 

AGd [D - l)rl -{D- 3)P 

In the Schwarzschild limit, i.e., I going to infinity, this formula turns into: 

{D - 2)/ii3_2rg-2 



4G 



< 



D 



(52) 



(53) 



(54) 



(55) 



(56) 



which is the well-known result for standard Schwarzschild BH, i.e., C < 0. 

Provided that the condition {Rq + f{Ro)) < is accomplished, the f{R) general case lead, like in the EH case, to 
C > for rjf > r^Q (the large BH region) and C < for < r^o (the small BH region). For ^ t^q (T close to 
To) C is divergent. Notice that in EH gravity, C < necessarily implies F > since Tq < Ti. 

In summary, f{R) theories accomplishing the condition: 



Ro + f{Ro) < 



(57) 
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provide a scenario analogous to the one described in full detail by Hawking and Page in fTD] for the EH case. This 
analysis was originally shown in [1121 1113j and then widely used in literature. 

In principle, particular f{R) models violating the condition Rq + /(i?o) < can be naively considered. However, in 
such cases, the mass and the entropy would be negative and therefore the AdS BH solutions for these models would 
be unphysical. Therefore Rq + f{Ro) < can be regarded as a necessary condition for f{R) theories in order to 
support the usual AdS BH solutions. Using Equation (27), this condition is fully equivalent to 1 + /'(i?o) > 0. This 
last condition has a clear physical interpretation in f{R) gravity theories (see |114| and references therein). Indeed, 
it can be interpreted as the condition for the effective Newton constant Ge// = Gd/(1 + /'(i?o)) to be positive. It 
can also be interpreted from the quantum point of view as the condition which prevents the graviton from becoming 
a ghost. Therefore the requirement Equation (57) can be also regarded as a thermodynamical viability condition for 
f{R) gravity theories when AdS configurations are studied [1121 1113] . 

Once that the main thermodynamical quantities have been derived, let us now revise the local and global stabilities 
phenomenology: For T < Tg, the only possible state of thermal equilibrium in an AdS space is pure radiation with 
negative free energy and hence there is no stable BH solutions. For T > Tq there exist two possible BH solutions; the 
small (and light) BH and the large (heavy) BH. The small one has negative heat capacity and positive free energy 
as the standard Schwarzschild BH. Therefore it is unstable under Hawking radiation decay. For the large BH there 
appear two possibilities; if Tq < T < Ti then both the heat capacity and the free energy are positive and the BH will 
decay by tunneling into radiation, but if T > Ti then the heat capacity is still positive but the free energy becomes 
negative. In this case the free energy of the heavy BH will be less than that of pure radiation. Then pure radiation 
will tend to tunnel or to collapse to the BH configuration in equilibrium with thermal radiation. This phenomenology 



will be applied to different f{R) models in Section VI 



B. BH Thermodynamics for KN Configuration 



Following the reasoning at the beginning of this section, the research concerning the KN configuration 
will be also done for BH with a well-defined horizon structure for negative values of Rq. In order to 
study the different thermodynamical properties of Kerr-Newman BH in f{R) theories, the usual approach 
is to start studying the temperature of the exterior horizon r^xt = ?'ext {Rq, a, Q, M). For that pur- 
pose, the Euclidean action method |122H125] is again considered. The procedure consists now of perform- 
ing the change of variables t — > —it, a — > za on the metric (38) in order to obtain the Euclidean section 



(see |120l 1121] for details). After some calculations, the inverse Hawking temperature now becomes: 



/3kn = 



47r (r 



£;,KN 



1 - 



Rq 
12 



(58) 



For null spin and charge, this expression becomes naturally Equation (45) 



BH horizon temperature can also be obtained through Killing vectors, as is explained in |126) where temperature 
is defined as: T^^kn = k/Att, with k the surface gravity defined now as ^fiXi' = '^Xv- It can be verified that k is 
the same at any horizon point and consequently T^^kn = Te,kn as obtained in |129j . 

Now that we know the expression for the temperature for KN configurations, the natural step in order to obtain 
the remaining thermodynamical quantities is to consider the corresponding extension of the Euclidean action with 
respect to the previously considered Equation (48). It now becomes: 



AS- 



E,KN — 



^ j d'x^\ [Ro + f{Ro) - F^^Fn 



(59) 



with y the integration region. As is described in 701, ^'^cl fully derived in |1201 1121] the result for IS.Se,kn becomes: 



BKN 



fi (Rq + f{Ro)) 

24 S 



12 a- 



2 1 



12 

Text 

Ra ) Rq Text 



^ ' ext 



(60) 



where is the electric potential of the horizon as seen from infinity: 



14 



and Q is the physical electric charge of the BH, obtained by integrating the flux of the electromagnetic fleld tensor 
at infinity, which happens to be: 



(62) 



From expressions Equations (61 1 and (62), it can be seen that $e and Q definitions are /(_R)-model independent. 
The reason for this lies in the fact that these calculations involve the vector potential and the electromagnetic field 
tensor, which are independent of the considered f{R) model. Further analysis of the action revealed jl20[ 1121] that 
it became singular for = 0, as could be expected from extremal BH, whose temperature kn = makes the 
Pkn factor diverge. Since thermodynamical potentials are obtained by dividing the action by the /3 factor, they still 
remain well-defined at Te,kn — 0. 

By using the expression Equation (60 1 we can immediately obtain Helmholtz free energy F. For rotating configu- 
rations, the correct definition is now: 

ASe 



F = 



VLhJ 



(63) 



where the term Q.^ J comes from the required Legendre transformation to fix angular momentum, being J the angular 
momentum of the BH and Q,h the angular velocity of the rotating horizon. For the explicit expressions, we refer the 
reader to [1201 1121] . Thus, performing some calculations one gets: 



F = (l + /'(i?o)) 



[36Q2 + 12rL+r,4.*i?o + a' (36- 



Ro)] 



24 Trrrt S 



(64) 



Once again, expression Equation (64) reduces to Equation (49) for a non-rotating uncharged BH. Provided that 



the condition 1 -(- f'{Ro) > is required to hold in order to obtain positive values of the mass, the inspection of the 

„ lim'i 
ext 



numerator of F renders that > for values of the horizon below r jt™'* (with an associated mass M through 



Equation (41 )), and < for larger values. 

Using the appropriate thermodynamical relations [I30j . we can derive the entropy S of the KN BH, which reads: 



S = l3{M-nHJ)-ASE = [1 + f'iRo)) 
If the area Ah of the exterior horizon is now computed, one gets: 



a') 



Ah = 47r {r 



(65) 



(66) 



Thus, if expression Equation (66) is substituted into the entropy expression Equation (65), the entropy can be 
expressed as: 

Ah 

4 



^ = (1 + f'iRo)) 



(67) 



Consequently 1 + /'(i?o) > is again a mandatory condition to obtain a positive entropy, as was shown in the AdS as 
well. Note that if for example we want to recover GR with a cosmological constant taking /(i?o) = ^2A, this relation 
turns into: S = Ah/4:, i.e., the famous Bekenstein result [128] . 

Once that the temperature T and the entropy S of the KN BH are known, the natural step further can be taken 
and study the heat capacity at constant scalar curvature Rq and at fixed spin a and charge Q parameters. From the 
definition: 



C = T 



df 



Ra.a.Q 



we obtain the expression: 



[a' (12 + rl,^ Ro) + 3 (4 4 rl 



■Ro)] 



Den{Q, r^xt, Ro, a) 



(68) 



(69) 



C - (l + /'(i?o)) 

with Den{Q, r^xt, Roi o.) = F)en defined as: 

Den = S [-36 rl,, + (-12 + r^* i?o) + 3 rt,, (4 + r^* i?o) - (3 + 12 r^* - 2 i?o)] 

At this stage, it is interesting to find out for which values of Rq, a, Q and M the denominator of the thermal 
capacity becomes zero, i.e., the thermal capacity goes through an infinite discontinuity, i.e., the BH experiences a 
phase transition. Following the analysis presented in |1201 1121] . two kinds of BH on this subject can be distinguished 
depending on the values of the a, Q and M parameters and scalar curvature Rq: 
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• fast BH, without phase transitions and always positive heat capacity C > 0. 

• slow BH, presents two phase transitions for two determined values of r^xt- 



In Figure [3] we have visualized the behavior of the temperature T and the heat capacity C of a BH for different 
values of mass parameter M , with fixed a, Q and Rq values. 
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Figure 3: For Ro = —0.2, we graphically display temperature (left) and heat capacity (right) of a BH as functions of the 
mass parameter M for the cases: (I) a — 0.5 y Q = : "slow" BH that shows a local maximum temperature Tmax and a 
local minimum temperature Tmin at the points where the heat capacity diverges, taking the latter negative values between 
Tmax y Tmin', (II) 1 ~ 0.965 y Q = : limit case where Tmax and Tmin merge, hence resulting in an inflection point in the 
temperature and an always positive heat capacity; (III) a = 1.5 y Q = : "fast" BH with both temperature and heat capacity 
monotonously growing (always positive too). 



As was remarked in |120l I12lj . KN-AdS BH, unlike Schwarzschild-AdS BH case [70], are allowed for any value of 
the temperature T. Hence stability of each BH will be exclusively given by the corresponding values of heat capacity 
C and free energy F as functions of the a, Q, M and Rq parameters, that ultimately define the BH. However, for 
a set of fixed values of a, Q and Rq, the mass parameter must be bigger than a minimum M""" (characterized by 
T = 0) to host BH configuration, otherwise radiation is the only possible equilibrium state up to such a minimum 
mass. For bigger masses, we shall distinguish between the fast and the slow BH. 

• Fast BH, with bigger values of the spin and the electric charge than the slow ones, shows a heat capacity 
always positive and a positive free energy up to a M = M value, and negative onwards. Thus, this BH is 
unstable against tunneling decay into radiation for mass parameter values of M < M For M > M 
free energy becomes negative, therefore smaller than that of pure radiation, that will tend to collapse to the BH 
configuration in equilibrium with thermal radiation. 

• Slow BH shows a more complex thermodynamics, being necessary to distinguish between four regions delimited 
by the mass parameter values: M < < AI^^ < M as follows: 

— For M™™ < M < and for < M < M '*™'*, both the heat capacity and the free energy are positive, 
which means that the BH is unstable and decays into radiation by tunneling. 

— If < M < M^^, the heat capacity becomes negative but free energy remains positive, being therefore 
unstable and decays into pure thermal radiation or to larger values of mass. 

— Finally, for M > M the heat capacity is positive whereas the free energy is now negative, thus tending 
pure radiation to tunnel to the BH configuration in equilibrium with thermal radiation. 

Let us conclude this section as we did the previous one, by emphasizing that, although not quantitatively, the 
thermodynamical behavior of KN~AdS BH in f{R) modified gravity theories is qualitatively analogous to that of GR 
[130) . This was an original important result presented for the first time in [1201 1121] . In that investigation, authors 
proved that some thermodynamical quantities like the physical mass M, the angular momentum J and the entropy 
S differ from those of GR by a (1 + /'(i?o)) factor. One might think that the free energy F and the heat capacity C 
also differ in that same multiplicative way, but because of the dependence with the normalized charge parameter Q, 
this same term (1 + /'(i?o)) leads to a more complex deviation from the GR behavior. 
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VI. THERMODYNAMICS IN AdS AND KERR-NEWMAN: PARTICULAR EXAMPLES 

A. Thermodynamics in AdS 



In order to illustrate the calculations sketched in Section V A let us consider some particular /(i?) models in whose 



heat capacity C and the free energy F will be determined that were originally studied in [1121 1113| . As was pointed 
out above, those two quantities provide the local and global stability of AdS BH. For the particular models that are 
studied below, the constant scalar curvature Rq could be calculated exactly by using Equation (27). For the sake of 
simplicity the D-dimensional Schwarzschild radius in Equation (31 1 was taken as Rs~^ = 2. The models we have 
considered are the following. 



Model I: f{R) = ai-RY 



This model has been widely studied because the aR^ term with a > can account for the accelerated expansion of 
the Universe |131| . This model can also explain the observed temperature anisotropics observed in the CMB, and can 
become a viable alternative to scalar field inflationary models, since reheating after inflation would have its origin on 
the production of particles during the oscillation phase of the Ricci scalar |132j . For (3 = 2, this model was recently 
pointed out [El [18] as a possible explanation to dark matter. 



Substituting this model in Equation (27 1 for arbitrary dimension and by only considering non- vanishing curvature 
solutions, the found solution is: 



2~D 



(70) 



For D > 2, the condition (2/3 — D)a < provides well-defined scalar curvatures Rq. Two separated regions have thus 
to be studied: Region 1 {a < 0, /3 > D/2} and region 2 {a > 0, /3 < D/2}. For this model we also get: 



1 + /'(i?o) = 



213- D 



(71) 



Notice that in region 1, 1 + f'{Ro) > for _D > 2, since in this case /3 > 1 is straightforwardly accomplished. In 
region 2, we find that for D > 2, the requirement Rq + f{Ro) < 0, i.e., 1 + /'(i?o) > 0, fixes f3 < 1, since this is the 
most stringent constraint over the parameter (3 in this region. Therefore the physical space of parameters in region 2 
is restricted to he {a > 0, (3 < 1}. 

In Figures |4]and[5j taken from [TT^ITTB] . one can see the physical regions in the parameter space (a, j3) corresponding 
to the different signs of (C, F). 



Model IV D=4, Regions 1&2 Model IV D=5, Regions 1&2 
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Figure 4: Thermodynamical regions in the (jaj, plane for Model I in _D = 4 (left) and D = 5 (right). 
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Mod el IV D = 10, Regio ns 1&2 
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Figure 5: Thermodynamical regions in the (|q|, |/3|) plane for Model I in D = 10. 



2. . Model II: f{R) = 



'R 

K I — 

a 



1 + /? 



R 



This model has been proposed in |106) as cosmologically viable and has attracted a lot of attention in order to 
constrain its validity in different astrophysical and cosmological scenarios. Throughout this study, let us consider 
n = 1 for the sake of simplicity. A vanishing curvature solution also appears in this model and the two nontrivial 
curvature solutions are given by |112[ I113j : 



[(k -2)D + A± y^VnD^ - 8D + 16 



2P{D~2) 



(72) 



where k > and k > (8-D — 1Q)/D'^ are required for real i?o solutions. The corresponding 1 + f'{Ro) values for 
Equation (72) are: 



l + /'(i?±) = 1 



4(D 



(VkD2 - SL* + 16 ± VkD) 



(73) 



Since 1 + /'(i?o) > is required, one has sign(i?Q ) = sign(a/3). In fact it turns out that 1 + /'(i?o ) is not positive for 
any allowed value of k and therefore this curvature solution Rq is excluded for our study (see [112( IllSj for details) . 
On the other hand, 1 + /'{Rq) > only requires k > for dimension D > 4 and therefore f'{Ro) < 0. Therefore 
only two accessible regions need to be studied: region 1 {a > 0, /3 < 0} and region 2, {a < 0, /3 > 0}. 



B. KN-AdS Thermodynamics 

In this section we will study the thermodynamics of the two f{R) models presented in the last section but when a 
Kerr-Newman configuration is considered. 

First, as was shown above, the condition 1 + /'(i?o) > must be satisfied in order to obtained physical results, so 
first of all the viability of each model depending upon the values of the parameters that define them will be checked 
out. Thus, for each model, we will study the range of parameters that allows the existence of Kerr-Newman BH, 
i.e., we will graphically display amax and amin (if present) in terms of those parameters, and study the region of a 
confined between these two surfaces. This phenomenology was absent in the SAdS study since the configuration was 
spherically symmetric and therefore a = 0. Finally we will focus on the thermodynamical quantities that define BH 
stability depending again upon the model range of parameters. In order to sketch graphically the thermodynamical 
behavior, let us then consider the BH parameter values: M = 1, a = 0.4 and Q = 0.2. Let us remind at this stage 
that this analysis is restricted to i?o < condition, as explained in the beginning of Section [V] 
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For the sake of simplicity, let us introduce the dimensionless variables notation [1201 1121) that we will use in the 
graphical visualization of the results to be presented: 



M 



a 

M 



a, ^^Q, RoM' 



Rq 



(74) 



where the previous symbols keep their original meaning as defined before. The considered models are as follows. 



Model I: f{R) = a{~Rf 



Let us remind the solution Equation (70) found above for this model. The viability condition 
1 + /'(i?o) > restricts the range of parameters that define this f{R) model. 

In Figure [6j we show the range of the spin parameter a for which BH are allowed, depending on the parameters a 
and /3 and for charge parameters values Q = and Q = 0.75. For the solution i?o < there exist two possible regions: 

• Region 3 {a < 0, /3 > 2}. For this region, the value of amax decreases suddenly from its normal value to near 
the frontier of the region, i.e., a « and /3 « 2; other values of a and /3 display a relatively low curvature and 
the existence of BH is assured. 

• Region 4 {a > 0, < 1}. This region only shows problems when < /3 < 1, where a„iax 0, but, as /3 becomes 
more negative, the surface of amax slowly acquires higher values, recovering its usual value for /3 = —2. 






Model I, Region 4, Q = 0' 




Figure 6: Model I: Region 3: {a < 0, /3 > 2}, and Region 4: {a > 0, < 1}. BH with a well defined horizon structure will 
only exist if they have a spin parameter below the upper surface amax. For the presented regions 3 and 4, the surface a-min 
does not exist. 

We have graphically schematized in Figure [7] the possible thermodynamical configurations of a BH for Regions 3 
and 4, i.e., Rq < following the results presented in jl201ll21] . As we can see in Region 3, only for values close to the 
region boundary (a « and /3 « 2) is BH stable both locally and globally (C > and F < 0). Also, as the absolute 
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value of a increases, it is more probable that the BH shows global instability (C > and F > 0) and finally global 
and local instability (C < and F > 0). In Region 4, however, we find the opposite behavior: near the limit of the 
region values (a « and /? « 1) we have global and local instability, and, as a gets bigger, the BH becomes more 
stable both locally and globally. 



Model I, Region 3 



Model I, Region 4 





-1.0 -0.5 0.0 0.5 

Log (or) 



Figure 7: Thermodynamical regions for Model I. We distinguish between three different regions: (i) C < and F > 0, in black; 
(ii) C > and F > 0, in gray; (iii) C > and F < 0, in white. 



2. . Model II: f{R) = -a 



'R 

K I — 

a 



1 + 



R 



Let us consider for the sake of simplicity again n = 1, thus having a two-parameters model, since we can define 
J = /3/a and then obtain: 



fiR) 



1 + -/R 



(75) 



Replacing the latter in Equation (27 1, one gets Equation (72), that for D = 4, and focusing on the negative constant 
curvature solution, yields: 



Fo 



1 



-K (1 — k) 



(76) 



Keeping in mind that we have to satisfy 1 + f'{Ro) > 0, k happens to be restricted to values of k > 1. On the 
other hand, computation of 1 + /'{Rq) reveals that Rq is only a valid solution for values of k and 7 in the Region 2: 
{k > 1, 7 < 0}, being as mentioned Rq < 0. In Figure [S] we show the range of the spin parameter a for which BH are 
allowed, depending on the parameters k and 7 for certain values of the charge Q. Region 2 displays an Umax surface 
that descends slightly as the value of k increases, and more significantly when 7 « 0, becoming when this limit is 
reached. 

We graphically schematize in Figure [9] the different possible thermodynamical configurations for the region in which 
Rq < 0, in this case the region 2. In this parameter space region, a combination of high values of k and small absolute 
values of 7 favor BH local and global stability, and, conversely, a combination of small values of k and high absolute 
values of 7 lead to BH instability, first globally and then also locally. 
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Figure 8: Model II. Region 1: {k > 1, 7 > 0}, and Region 2: {k > 1, 7 < 0}. BH with a well-defined horizon structure will 
only exist if they have a spin parameter below the upper surface amax, and above a second surface amin (in case it exists, only 
in region 1 for this model) for certain values of k and 7. 



VII. COSMOLOGICAL SOLUTIONS IN MODIFIED GRAVITY 



Modified gravity has become very popular in recent years because it is capable to reproduce the dark energy epoch, 
and even the inflationary phase, with no need to introduce an additional component in the universe. Hence, the 
reconstruction of appropriate theories that contain some particular cosmological solutions have involved an important 
effort (see [32l|38l[42-46, 101t .l33. - .149j ). The same is true for the reconstruction of so-called viable modified gravities, 
which can recover GR in some limit while they introduce additional effects at cosmological scales (see |106H110) ). 
However, as in every dark energy model, the required fine tuning on the form of the action as on free parameters of 
the model have resulted in a large effort in searching additional and independent predictions in the frame of modified 
gravities, including other aspects different than cosmology as the ones studied in the above sections (see also the 
reviews |26H33j and references therein). In this section we are interested to make a short review on cosmological 
solutions in the frame of modified gravities, and in particular in F{R) gravity, while the thermodynamical aspects are 
studied in the following sections. As it is widely accepted, our universe seems to be well described by FLRW metrics: 



1 

1 — kr^ 



r^de^ + sin^ 9 dtp' 



(77) 



and specifically the universe seems to be spatially fiat, fc = 0. Hence, by introducing the metric in F{R) field 
Equation ([s]), modified FLRW equations are obtained: 



-3H^ -2H = 



Fn 



3Fr 



RFr - F 



— SHRFffji 



k'p^ + R^Frrr + 2HRFrr + RFrr + -{F- RFr) 



(78) 
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Figure 9: Thermodynamical regions with negative scalar curvature i?o < of model I. We distinguish between three diflterent 
regions: (i) C < and F > 0, in black; (ii) C > and F > 0, in gray; (iii) C > and F < 0, in white. 



Recent astronomical data indicates that the well-known ACDM model where DE is represented by an effective cosmo- 
logical constant fits quite well the observational data. Nevertheless, it is shown that F{R) gravity can even reproduce 
exactly the ACDM model in absence of a cosmological constant s shown below (see [571ll50| ). Let us consider ACDM 
cosmology: 



K 



ff2 I ^ „-3„-3Af 



(79) 



where we have introduced the number of e-foldings N = In-^ instead of the cosmological time t, and we have assumed 



a pressureless fluid. We can easily solve the continuity equation for pm'- 



(80) 



while the first FLRW equation in Equation (78 1 can be rewritten in terms of the number of e-foldings, which yields 



= -9G {N (i?)) (4G" (N (R)) + G" {N (R))) + (^3G (N (R)) + ^G' (N {R))^ 



(81) 



where G{N) — H^. Hence, for ACDM cosmology Equation (79), the scalar curvature is given by i? = 3G'{N) 
12G{N) = 12i/Q + poaQ^e""^^ , and the Equation (81 1 gives a hypergeometric equation [HTj: 



= 3(i?-9i7^)(i?-12i/„^)^^ 



(82) 



whose solution is given by: 

F{x) = AF{a, (3, 7; x) + Bx^-^F{a - 7 + 1, /? - 7 + 1, 2 - 7; x) 
where x — — 3, {A, B} are integration constants, and the parameters are: 



(83) 



7 = + = , a/3 : 



(84) 
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Hence, the action Equation ( 83 ) is capable to reproduce ACDM model with no cosmological constant. Nevertheless, the 



action Equation (831 is composed by hypergeometric functions, which do not seem to provide a successful explanation 



for dark energy, since the theory turns out much complicated than the usual GR with a cosmological constant. 
However, other viable models of F{R) gravity capable to mimic a cosmological constant at the present time have been 
suggested in the literature (see |106H110] ). which in general seem to cross the phantom barrier |151| . Then, let us now 
study a model where the dominant component is phantom-like. Such kind of system can be easily expressed in the 
standard General Relativity when a phantom fluid is considered, where the FLRW equation reads (t) — ^ pph (the 
subscript ph denotes the phantom nature of the fluid). As the EoS for the fluid is given by Pph — WphPph with < — 1, 
by using the conservation equation, the solution yields a{t) — ao(ts — ^ where aq is a constant, — — ^i^^ ^ ) 
and tg is the so-called Rip time. Then, the solution describes a universe that ends in a Big Rip singularity at the time 
tg. The same behavior can be achieved in F{R) theory with no need to introduce a phantom fluid. The expression for 
the Hubble parameter as a function of the number of e-foldings is given by H^{N) = Hle^^'"°. Then, the Equation 
(81), with no matter contribution, takes the form: 



AR^EM+BF{R)=0 



dR^ ' dR 

where A = -Ha{l + Hq) and B = il±Mo) . This equation is an Euler equation whose solution yields: 



F{R) = Cii?'"+ + C2R" 



where m± 



1- A±y/{A-l)^ -AB 



(85) 



(86) 



Thus, the phantom cosmology a{t) — ao(is — can be also obtained in the frame of F{R) theory and no phantom 
fluid is needed. 

In addition, the equivalence between f{R) gravity and scalar-tensor theories Equation ^ can be used to easily 
obtain cosmological solutions in f{R) gravity. For the metric Equation (77), Friedmann equations read: 



3H^^+3ij2p(0) + iQ(0)-p„ 







d^Pjcl)) 
dt^ 



(87) 



We redefine the scalar field such that it is chosen to be the time coordinate cj) = t. Hence, taking into account the 



Equations (87) and the continuity equation for the perfect fluid pm, the Hubble parameter may be calculated as a 



fimction of the scalar field cj), H = g{4>). By combining the Equation (87), the function Q{4>) is deleted, and it yields 



d'^PU) , ..dPidi] 

- 9^'^) + 2g'(0)^('^) + (1 + w„Op,„o exp 



-3(l + w;„) / d^g{^) 







By solving this equation for a given function P{(t>), a cosmological solution H{t) is found, and the function (3(0) is 
obtained by means of equations in Equation ( 87 ) : 



QW = -6(.g(0))2p(</>) - 6g(0) 



d(j) 



2p™o exp 



-3(1-1- w™) / d^gi^) 



(89) 



If we neglect the contribution of matter, then the Equation (88) is a first order differential equation on gi^cj)), and it 



can be easily solved. The solution found is the following |147 



g(0) = C/pM-/pM 



2F3/2(0) 



(90) 



where C is an integration constant. Hence, we have shown that F{R) gravity represents a serious candidate to explain 
the phenomena of late-time acceleration in the expansion of the universe. In the following sections, we explore the 
thermodynamical properties of FLRW metrics in the frame of modified gravity as well as an attempt to generalize 
Cardy-Verlinde formula. 



VIII. FIRST LAW OF THERMODYNAMICS AND FLRW EQUATIONS 



In this section, we review some of the most important results obtained in the study of FLRW metrics and their 
thermodynamics, specifically the relation between the first law of thermodynamics and FLRW equations, which can 
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be deduced by assuming the well known relation between the entropy and the horizon area, as well as the relation 
between the gravity surface and the temperature [82H84] . Such study has been extended not only to f{R) gravity, 
as is also shown below (see [Hi]), but also to many other different scenarios In fact, it was found in [SO] 

and extended in [811 for f{R) gravities, a direct deduction of gravity field equations by assuming the relation of the 
horizon area and the entropy. Nevertheless, here we are only interested to show the deduction of FLRW equation 
from the thermodynamics defined in the corresponding horizons of a particular spacetime. Let us consider a ?i + 1 
FLRW metric Equation (771, with arbitrary spatial curvature fc = — 1, 0, 1. This metric may be rewritten as: 



ds^ = habdx'^dx'' + f^dnl_^ 



(91) 



where the indexes a, 6 = 0, 1, being a;° = t, = r, and f — a{t)r, while drt^_i expresses the line element of an n — 1 
unit sphere. Then, the apparent horizon for this general kind of spacetimes is defined as h'^'^dafdbr = 0, whose radius 
yields: 



1 



(92) 



As will be shown below, one can assign an entropy proportional to the area of the apparent horizon as well as an 
energy flux through the horizon. This is the reason for considering the apparent horizon instead of other horizons 
defined in FLRW spacetimes. We may define the work density and the energy flux vector as (82n84) : 



w 



*a = T%r + Wdaf 



(93) 
(94) 



where the work w refers to the change experimented by the apparent horizon, while the vector refers to the flux 
of energy through the horizon. Then, the variation of the internal energy can be written as: 



S/E = A^ + 'wS/V 



(95) 



where A is the area and V is the volume of an n-dimcnsional space with radius f. This equation represents the unified 
first law according to [152J. The first term in the r.h.s. in Equation (95) takes the form: 



8ttG 



E 



where k is the surface gravity: 



Note that the energy is defined as: 



-.da (^V^h'^'dbf) 



167rG 



(1 - h-'dafdbf) 



(96) 



(97) 



(98) 



Hence at the apparent horizon, the second term in the r.h.s. 



horizon, so that the expression Equation ( 96 1 can be identified with the thermodynamical relation dE 



in Equation (96) vanishes, as on the black hole event 

-TdS if 



the surface gravity is assumed to be the temperature k/2t:, and the entropy yields the usual expression, proportional 
to A/ AG. Now if we assume a perfect fluid filling the universe, with an energy-momentum tensor given by T^i/ — 
u, the energy fiux vector yields: 



Then, the change of the internal energy through the apparent horizon Equation ( 96 ) becomes: 

dE = A^ = -A{p + p)HrAdt 
Identifying the entropy and temperature associated to the apparent horizon, 

A „ 1 



S 



AG 



T 



2TTfA 



(99) 



(100) 



(101) 
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and using the first law of thermodynamics dE = —TdS, the second FLRW equation is obtained: 



SttG ^ . . . k 



-{p + p) + ^ (102) 



By combining it with the continuity Equation (80), the first FLRW equation is also obtained. Hence, FLRW equations 
can be deduced by using only thermodynamical arguments and identifying the entropy of the apparent horizon and 
its temperature, which is an identity that supports the idea of the relation between gravitational field equations and 
thermodynamical ones, as suggested by Jacobson in [50] . 

In the same way, this analysis can be easily extended to f{R) gravity as well as to other modified theories of gravity. 
Specifically, for F{R) gravity, the entropy relation is defined (see previous sections) as follows: 

AF'(R) , , 



Hence, following the same steps performed above for GR, the FLRW equations for f{R) gravity Equation (781 are 
recovered j85j . Then, a strong relation between thermodynamics of the horizons and gravitational field equations is 
established. In the next section, another aspect of the relation of FLRW equations and thermodynamics is explored, 
in that case following the holographic principle through the entropy bounds and the Cardy-Verlinde formula. 



IX. GENERALIZATION OF CARDY-VERLINDE FORMULA 



This section is devoted to consideration of the Cardy-Verlinde formula for general scenarios, and in par- 
ticular for modified gravity. The Cardy-Verlinde formula was established by Verlinde in |98j as a ther- 
modynamical relation between FLRW equations with a conformal fiuid (a radiation-like fluid) and the 
Cardy formula (see [57]), and it was extended for a general perfect fluid in [99. In |100j . attempts 
to reproduce the Cardy-Verlinde formula was suggested for different scenarios, including modified gravity. 
Other aspects of the Cardy-Verlinde formula such as the presence of a viscous fluid have been studied 
in |153H155] . In this section, we mainly review those results. We consider a (n ~{- l)-dimensional spacetime de- 



scribed by the FLRW metric Equation (77). By inserting the metric Equation ([77| in the Einstein field equations, 
the FLRW equations for GR are derived: 

n n — 1 — ' n — — ' 

^ ' i=i 



=1 



Here pi = Ei/V and pi are the energy-density and pressure of the matter component i that fills the Universe. In this 
section, we consider only the k = 1 closed Universe. Moreover, we assume an equation of state (EoS) of the form 
Pi = WiPi with Wi constant for each fluid, and assume no interaction between the different components. Then, the 



continuity Equation (80) has the form now pi + nH {pi + Pi) =0 and by solving it, we find that the i fluid depends 
on the scale factor as: 

p, cx (105) 

Let us now review the case of |99j . where just one fiuid with EoS p = wp and w = constant is considered. The total 
energy inside the comoving volume V , E = pV, can be written as the sum of an extensive part E-e and a subextensive 
part Eq, called the Casimir energy (for other ways to account the Casimir energy in a FLRW universe, see [15611157] ). 
and takes the form: 

E{S,V)^Ee{S,V) + ^Ec{S,V) (106) 

Under a rescaling of the entropy {S — > XS) and the volume {V — > \V), the extensive and subextensive parts of the 
total energy transform as: 

^e(A5, XV) = XEe{S, V) , EciXS, XV) = X^-^/^'EdS, V) (107) 

Hence, by assuming that the Universe satisfies the first law of thermodynamics, the term corresponding to the Casimir 
energy Eq can be seen as a violation of the Euler identity according to the definition in ^98j : 



Ec = n{E+pV -TS) 



(108) 
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Since the total energy behaves as ~ a""'^ and by the definition Equation (106), the Casimir energy also goes 
as Ec ^ a""™. The FLRW universe expands adiabatically, dS = 0, so the products Eqci^^ and i?Eo"™ sho uld b e 
independent of the volume V, and be just a function of the entropy. Then, by the rescaling properties Equation ( 107 1, 
the extensive and subextensive part of the total energy can be written as functions of the entropy only |99| : 

47ra"'" 27ra""' ^ ^ 



Here a and /3 are undetermined constants. By combining these expressions with Equation (106), the entropy of the 
Universe is written as a function of the total energy E and the Casimir energy Eq 

S = ^--^^Ec{2E-Ec)^ (110) 

which for 1^1 = 1/71 (radiation- like fluid) reduces to [55] : 

27ra 



/a/3 



^Ec{2E-Ec) (111) 



which has the same form as the Cardy formula given in [57] • The first FLRW Equation (104 1 can be rewritten as a 
relation between thermodynamics variables, and yields: 

= —a^EBii{2E-EBH) , where 5H = (n-l)^, Ebh ^ n{n ~ 1) . (112) 



It is easy to check that for the bound proposed in [98], Eq < Ebk, the equation for the entropy Equation (111) 
coincides with the first FLRW Equation (112) when the bound is reached. We will see below that when there are 
several fiuid components, the same kind of expression as in |98| cannot be found as well as its extension to modified 
gravity, nor is there the same correspondence with the FLRW equation when the bound is saturated. 



A. Multicomponent Universe 

If m fluids are considered with arbitrary EoS, Pi — WiPi, the expression for the total entropy is simple to derive just 
by following the same method as above. The total entropy is given by the sum of the entropies for each fluid [100) : 

5 = f]5. = r^^^^E^cm - E^c)] (113) 

This expression cannot be reduced to one depending only on the total energy unless very special conditions on the 
nature of the fluids are assumed. Let us for simplicity assume that there are only two fluids with EoS given by 
Pi — wipi and p2 — W2P2, wi and W2 being constants. We can substitute the fluids by an effective fluid described by 
the EoS: 

P1+P2 W2-W1 

PcS = WoffPoff , where WcS = ; = wi + — — (114) 

P1 + P2 I + P1/P2 



and Pcff — -^{pi + P2), PeS = ^{pi + P2)- Then, by using the energy conservation Equation (80 1, we find pi ~ 
(a/ao)""(i+"'i) and p2 (a/ao)-"(i+'"=), where Oq is assumed to be the value of the scale factor at the time Iq. The 
effective EoS parameter Wcs can be expressed as a function of the scale factor a(t): 

WcS = Wl + —, r (115) 

The total energy inside a volume V becomes: 

Et = Ei+E2 0^ (a/ao)-'™i + (a/ao)~"™' (116) 

As the energy is proportional to two different powers of the scale factor a, it is not possible to write it as a function 
of the total entropy only. As a special case, if the EoS parameters are wi — W2 — w^g, the formula for the entropy 



reduces to Equation (110 1, and coincides with the CV formula when Wcs = 1/n. 
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As another case, one might consider that for some epoch of the cosmic history, wi ^ W 2 ■ T aking also a >> ao, 
we could then approximate the total energy by the function Et oc a~"™^ . From Equation ( 106 ) the Casimir energy 
would also depend on the same power of a, Eq oc a^'™^. The expression Equation (110 1 is again recovered with 

w = ^2- 

Thus in general, when a multicomponent FLRW Universe is assumed, the formula for the total entropy does not 
resemble the Cardy formula, nor does it correspond to the FLRW equation when the Casimir bound is reached. It 
becomes possible to reconstruct the formula Equation ( 111 ), and establish the correspondence with the Cardy formula, 
only if we make specific choices for the EoS of the fluids. 



B. Inhomogeneous EoS Fluid and Modified Gravity 



Let us now explore the case of an n + 1-dimensional Universe filled with a fluid satisfying an inhomogeneous EoS. 
This kind of EoS, generalizing the perfect fluid model, has been considered in several papers as a way to describe 
effectively the DE (see |158L 1159] ). We assume an EoS expressed as a function of the scale factor: 



p = w{a)p + g{a) 



(117) 



This EoS fluid could be taken to correspond to modified gravity, or to bulk viscosity [1581 1159j . By introducing 
Equation ill7l in the energy conservation Equation (|80| we obtain: 



, , , n(l + w(a)) , , 

p'{a) + ^ —p{a) 

a 



5(a) 



(118) 



Here we have performed a variable change t — t{a) such that the prime over p denotes derivative with respect to the 
scale factor a. The general solution of this equation is: 



p{a) 



c 



where F{a) 



l + w{a') 



da 



(119) 



and C is an integration constant. As shown above, only for some special choices of the functions w{a) and g{a) can 
the formula Equation (fill be recovered. Let us assume, as an example, that w(a) = —1 and g{a) = —a™, with m 
being constant. Then, the energy density behaves as p oc a™. Hence, by following the same steps as described above, 
the extensive and subextensive energy go as a™"*"", and by imposing conformal invariance and the rescaling properties 



Equation (107), we calculate the dependence on the entropy to be 



Ev. — 



Ec = 



^— {2n-\-m / n) 



(120) 



The expression for the entropy is easily constructed by combining these two expressions and substituting the extensive 
part by the total energy. This gives us the same expression as in Equation (110) with w = —{n + m)/n. Note that 
for m — — (1 + n), the formula Equation (111) is recovered and also its correspondence with the CFT formula. 
However for a generic power m, the CV formula cannot be reconstructed, like the cases studied above. Only for some 
special choices does the correspondence work, leading to the identification between the FLRW equation and the Cardy 
formula. 

In the same way, we may reconstruct Cardy- Verlinde formula for modified gravity, and specifically for f{R) gravity. 
Note that f{R) gravity can be seen as GR plus an inhomogeneous fiuid. For a closed 3 + 1 FLRW Universe, the 
modified FLRW equations are expressed now as: 



1 



F{R) -~ 3{H^ + H)F'{R) + 3HF"{R)R = p„ 



-F(R) 



H 



F'{R)~[{duF'{R))+2H{dtF'{R))] 



K^Pn 



(121) 



where primes denote derivatives respect to R and dots with respect to t. These equations can be rewritten in order 



to be comparable with those of standard GR, similarly as in Equation ( 78 ) , by shifting terms to the right side of 



Equation (121). For such a propose, the geometric terms can be presented as an effective energy-density Pf(r) £^nd a 
pressure PF{Ry- 



Pf(r) 



1 



Pf{r) 



1 



K^F'{R) 



R^F"'{R) 



K^F'{R) 
2HRF"{R) 



RF'{R) - F{R) 



- 3HRF"{R) 



i?F"(i?) + -{F{R) - RF'{R)) 



(122) 
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Then, an EoS for the geometric terms can be defined as Pf{r) = Wp(^ii)PF{B:) ■ We can define an effective energy-density 
P = Pm/F'{R) + pF{R) and pressure p = Pm/F'{R) +Pf{r)- Hence, for some special cases the formula for the entropy 
developed above can be obtained in F(i?)-gravity (for an early attempt deriving a CV formula in a specific version of 
F(i?)-gravity, see I 160j ). For example, for an F{R) whose solution gives p oc a-3(i+iuett) ^ |;jjg formula for the entropy 
(110 1 is recovered although in general, as in the cases studied above, no such expression can be given. On 
land, one could assume that the geometric terms do not contribute to the matter sector. Supposing a 
constant EoS matter fluid, the expression for the entropy is given by Equation (110), although the cosmic Cardy 



Equation 
the other 



formula Equation (112 1 does not have the same form, and in virtue of the modified first FLRW Equation (121) the 



form of the Hubble entropy ^h, the total energy E and the Bekenstein energy Ebh will be very different. It is not 
easy to establish correspondence between two such approaches. Note that using the effective fiuid representation the 
generalized CV formula may be constructed for any modified gravity. 

Let us now consider de Sitter space solution in F(i?)-gravity (for review of CV formula in dS or AdS spaces, see 
[1611 1162] ). As was pointed in [55], almost every function F{R) admits a de Sitter solution. This can be easily seen 
from the first FLRW equation in Equation (121). A de Sitter solution is given by a constant Hubble parameter 
H{t) = Hq; then by inserting in Equation (121 1 we obtain the following algebraic equation: 



FjRo) 
2F'(i?o) 



(123) 



Here Rq ~ 12i?Q and the contribution of matter is neglected. Then, for positive roots Hq of this equation, the 
corresponding F{R) leads to the de Sitter solution which may describe infiation or DE. In this case the formula for 
the entropy Equation ( |110 ) can be reproduced for w = — 1, and even the universal bound Equ ation (|125|) can hold by 



taking a critical size of the Universe. The formula that relates the cosmic bounds in Equation ( 128 1 is easily obtained 
also in F(i?)-gravity for a de Sitter solution. In such a case one can identify: 



5'h — 



2G 



aV F{R„) 
24G F'{Ro) 



Sbh — 



V 
2Cki 



- <?2 



(124) 



Thus, one can conclude that 



which corresponds to the first FLRW equation written as 5^ + {S^, — S-qb)' 
dynamical entropy bounds are not violated for modified gravity with de Sitter solutions. Note that quantum gravity 
effects may be presented also as an effective fluid contribution. In case when de Sitter space turns out to be the 
solution, even with the account of quantum gravity the above results indicate that dynamical cosmological/entropy 
bounds are valid. In other words, the argument indicates the universality of dynamical bounds. It seems that their 
violation is caused only by future singularities if they are not cured by quantum gravity effects. Note that a large 
number of modifled gravity theories do not contain future singularities; they are cured by higher derivatives terms. 



X. ON THE COSMOLOGICAL BOUNDS NEAR FUTURE SINGULARITIES 



In [98j . Verlinde proposed a new universal bound on cosmology based on a restriction of the Casimir energy Eq] 
cf. his entropy formula Equation (111). This new bound postulated was: 



Ec < E'en 



(125) 



where -Ebh = n(n — 1) ^^^i ■ It was deduced by the fact that in the limit when the Universe passes between 
strongly and weakly self-gravitating regimes, the Bekenstein entropy S-q = '^^F and the Bekenstein-Hawking entropy 
5'bh = {n — 1)4^, which deflne each regime, are equal. This bound could be interpreted to mean that the Casimir 
energy never becomes able to reach sufficient energy, Ebh, to form a black hole of the size of the Universe. It is easy 
to verify that the strong {Ha > 1) and weak {Ha < 1) self-gravity regimes have the following restrictions on the total 
energy: 



F < Ebh for Ha<l 
F > Ebh for Ha>l 



(126) 



From here it is easy to calculate the bounds on the entropy of the Universe in the case when the Verlinde formula 
Equation (111) is valid; this is (as shown in the above sections) for an effective radiation dominated Universe WcS ^ 
1/n. The boimds for the entropy deduced in for /c = 1 are: 



S < Sb for Ha < 1 
S > Sb for Ha > 1 



(127) 
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where Sb is the Bekenstein entropy defined above, and iSh is the Hubble entropy given by Equation (112 1. Note 



that for the strong self-gravity regime, Ha > 1, the energy range is Ec < Eqh < E. According to the formula 



Equation (111) the maximum entropy is reached when the bound is saturated, Eq = i?BH- Then S = Sn, such that 
the FLRW equation coincides with the CV formula, thus indicating a connection with CFT. For the weak regime. 
Ha < 1, the range of energies goes as Ec < E < Ebh and the maximum entropy is reached earlier, when Ec = E, 
yielding the result S = Sb- The entropy bounds can be extended to more general cases, corresponding to an arbitrary 



EoS parameter w. By taking the bound Equation (125) to be universally valid, one can easily deduce the new entropy 



bounds for each regime, from the expression of the entropy Equation (110). These new bounds, discussed in 



differ from the ones given in Equation (127), but still establish a bound on the entropy as long as the bound on 
Ec expressed in Equation ( 125 ) is taken to be valid. The entropy bounds can be related through the first FLRW 
equation, yielding the following quadratic expression (for k — 1): 

+ {Sb - SBnf = Si (128) 



Here, we would like to review what happens to the bounds, particularly to the fundamental bound Equation ( 125[ ), 



when the cosmic evolution is close to a future singularity (see |100j ): then the effective fluid dominating the cosmic 
evolution could have an unusual EoS. As shown below, for some class of future singularities, such a bound could 



soften the singularities in order to avoid violation of the universal bound Equation (125). It could be interpreted to 
mean that quantum effects become important when the bound is reached. However, as the violation of the bound 
could happen long before the singularity even in the presence of quantum effects, it could be a signal of breaking of 
the universality of the bound Equati on (|125[ ). Let us first of all give a list of the possible future cosmic singularities, 
which can be classified according to |163j as: 

• Type I ("Big Rip"): For t ^ ts, a — ?> oo and p — > oo, \p\ — > oo. 

• Type II ("Sudden"): For t ^ t^, a ^ a^ and p ^ ps, |p| — > cx). 

• Type III: For t ^ tg, a as and p — )• oo, \p\ — > oo. 



• Type IV: For t — !■ , a — > Os and p — > ps, p Ps but higher derivatives of Hubble parameter diverge (see |164) ). 

Note that the above list was suggested for the case of a flat FLRW Universe (see [165H169] ). As we consider in this 
section a closed Universe {k = 1), we should make an analysis to see if the list of singularities given above is also 
valid in this case. It is straightforward to see that all the singularities listed above can be reproduced for a particular 
choice of the effective EoS. To show how the cosmic bounds behave for each type of singularity, we could write an 
explicit solution of the FLRW equations, expressed as a function of time depending on free parameters that will be 
flxed for each kind of singularity. Then, the Hubble parameter may be written as follows: 



Hit)^J^^^p-\=H,its-ty^+Ho (129) 
y n(n — 1) 

where m is a constant properly chosen for each type of singularity. Note that this is just a solution that ends in 
the singularities mentioned above, but there are other solutions which also reproduce such singularities. We will 
study how the cosmic bounds behave near each singularity listed above. As pointed out in [163. .170 190 . around a 
singularity quantum effects could become important as the curvature of the Universe grows and diverges in some of 
the cases. In other words, approaching the finite-time future singularity the curvature grows and universe reminds 
the early universe where quantum gravity effects are dominant ones because of extreme conditions. Then one has 
to take into account the role of such quantum gravity effects which should define the behavior of the universe just 
before the singularity. Moreover, they may act so as to prevent the singularity occurrence. In a sense, one sees 
the return of quantum gravity era. However, the consistent quantum gravity theory does not exist so far. Then, in 
order to estimate the influence of quantum effects to universe near to singularity, one can use the effective action 
formulation. We will apply the effective action produced by conformal anomaly (equivalently, the effective fluid with 
pressure/energy-density corresponding to conformal anomaly ones) because of several reasons. It is known that at 
high energy region (large curvature) the conformal invariance is restored so one can neglect the masses. Moreover, one 
can use large N approximation to justify why large number of quantum flelds may be considered as effective quantum 
gravity. Finally, in the account of quantum effects via conformal anomaly we keep explicitly the graviton (spin 2) 
contribution. The conformal anomaly Ta has the following well-known form: 

TA^b(F+ ^Di? j + b'G + b"UR (130) 
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Here we assume for simplicity a 3 + 1 dimensional spacetime. Then, F is the square of a 4D Weyl tensor and G is 
the Gauss-Bonnet invariant: 



kl 



(131) 



The coefheients b and b' in Equation (130) are described by the number of N scalars, N1/2 spinors, A'^i vector fields, 
A'2 gravitons and A^hd higher derivative conformal scalars. They can be written as: 



N + + 12Ni + 611A^2 - SiVHD 

120(47r)2 



b' 



N + lliVi/2 + 62Ni + 14117V2 - 28A^HD 
360(47r)2 



(132) 



As b" is arbitrary it can be shifted by a finite renormalization of the local counterterm. The conformal anomaly Ta 
can be wr itten as Ta = —Ra + 3p^, where pA and pA are the energy and pressure densities respectively. By using 
Equation ( 130 1 and the energy conservation equation pA + 3H(pA +Pa) = 0, one obtains the following expression for 

PA mm^ 



PA 



1 



1 



dta^HTA 



dta^H 



-126ij2 + 246'(-ij2 + H'^H + H"^) 



-(46 + Qb"){H + IHH + AH^ + 12H^H) 
The quantum corrected FLRW equation is given by: 

i/2 



8ttG , , 1 

-^\P^ Pa) ^ 



(133) 



(134) 



We study now how the bounds behave around the singularity in the classical case when no quantum effects are added, 
and then include the conformal anomaly Equation ( 130 1 quantum effects in the FLRW equations. We will see that 
for some cases the violation of the cosmic bound can be avoided. 



A. Big Rip Singularity 



This type of singularity has been very well studied and has become very popular as it is a direct consequence in 
the majority of the cases when the effective EoS parameter is less than —1, the so-called phantom case |170H190] . 
Observations currently indicate that the phantom barrier could have already been or it will be crossed in the near 
future, so a lot of attention has been paid to this case. It can be characterized by the solution Equation (1291 with 
TO < — 1, and this yields the following dependence of the total energy density on the scale factor near the singularity, 
when a ^ 1, for a closed Universe (fc = 1): 



n(n — 1) 
IGTrG 



if 2 



1 



-n(l+to) 



for t^i. 



(135) 



where we have chosen iJi = 2/n|l + w| with w<— 1,to = — 1 and iJo = for clarity. This solution drives the 
Universe to a Big Rip singularity for t ^ t^, where the scale factor diverges. If the singularity takes place, the bound 
Equation (1251 has to be violated before this happens. This can be seen from Equation (1351, as the Casimir energy 
behaves as Eq oc a"''"' while the Bekenstein-Hawking energy goes as Ebh oc a"~^. Then, as w < —1, the Casimir 
energy grows faster than the BH energy. Therefore, close to the singularity where the scale factor becomes very big. 



the value of Eq will be much larger than Eq, thus violating the bound Equation (125). Following the postulate 
from fSH] one could interpret the bound Equation ( 125 ) as the limit where General Relativity and Quantum Field 
Theory converge, such that when the bound is saturated quantum gravity effects should become important. QG 
corrections could help to avoid the violation of the bound and may be the Big Rip singularity occurrence. As this is 
just a postulate based on the CV formula, which is only valid for special cases as shown in the sections above, the 
bound on Eq could not be valid for any kind of fluid. 



Let us now include the conformal anomaly Equation (130) as a quantum effect that becomes important around 



the Big Rip. In such a case there is a phase transition and the Hubble evolution will be given by the solution of the 
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FLRW Equation (134). Let us approximate to get some qualitative results, assuming 3 + 1 di mensions. Around tg 
the curvature is large, and \pa\ >> {3/k^)H^ + k/a?. Then p ~ —pA^ and from Equation (133)we get: 

p + AHp = H [-12&ij2 + 246'(-ij2 + h'^H + H^) - (46 + Qh"){H + 7HH + 4ij2 + 12H^H) 

We assume that the energy density, which diverges in the classical case, behaves now as: 



(136) 



(137) 



where A is some negative number. By using the energy conservation equation p + 3-ff(l + w)p = 0, the Hubble 
parameter goes as H ^ l/{ts — t). We can check if this assumption is correct in the presence of quantum effects by 
inserting both results in Equation (136). We get: 



3iJ''(-13fe + 24fe') 



(138) 



Hence as 6 > and b' < 0, p becomes negative, which is an unphysical result. Thus p should not go to infinity in 
the presence of the quantum correction. This is the same result as obtained in |163) where numerical analysis showed 
that the singularity is moderated by the conformal anomaly, so that the violation of the bound that naturally occurs 
in the classical case can be avoided/postponed when quantum effects are included. 
Now we consider the case where f{R) has the form: 



(139) 



when the curvature is small or large. Then if the matter has an EoS parameter w > —1, by solving ( 121 1 we find: 

t^(^+^) when 



n(l+if;) 



> 



Then there may appear a singularity at t 



(i,-t)'.a+») when^^<0 
0, which corresponds to the Big Bang singularity, or at t 



(140) 



corresponds to the Big Rip singularity. Since the Casimir energy behaves as Eq 



but Ebh 



= tg, which 
^, only when 



t — >■ 0, Eq dominates in case that n > 2 and w > or in case that n > 2 and w > 0. Even in the phantom phase 
where 



2a 



n{l-\-w) 



< 0, the bound Equation ( 125 ) is not violated. 



B. Sudden Singularity 



This kind of singularity is also problematic with respect to the bounds. Nevertheless, as long as the energy density 
p does not diverge, the violation of the bound may be avoided for some special choices. The sudden singularity can 
be described by the solution Equation (129) with < m < 1, and constants Hq i > 0. Then the scale factor goes as: 



a{t) oc exp 



Hi 
m + 1 



its - t)"+' + Hot 



(141) 



which gives a{t) ~ e (de Sitter) close to ts- From the first FLRW equation the total energy density becomes: 



p = H\t) + — - [Hi{ts ~ t)™ + HoY + exp 



2^%-tr+' 

m + 1 



2HQt 



which tends to a constant p • 



Hi 



for t ts- Then the Casimir energy grows as Eq oc iJpa" 



a" 



(142) 



while 



ii^BH oc a"~ close to ts- The BH energy grows slower than the Casimir energy, and the bound is violated for a finite 
t. However, by a specific choice of the coefficients, the violation of the bound Equation ( 125 ) could be avoided. For 
Hq = 0, and by some specific coefficients, the bound could be obeyed. In general, it is very possible that Ec exceeds 
its bound. In the presence of quantum corrections, the singularity can be avoided but the bound can still be violated, 
depending on the free parameters for each model. We may assume that in the presence of the conformal anomaly for 
n = 3, the energy density grows as |163j : 



p^ po + pi{ts - t) 



X 



(143) 



where po and pi are constants, and A is now a positive number. Then the divergences on the higher derivatives of the 
Hubble parameter can be avoided, as is shown in |163j . Nevertheless, Ec still grow s fas ter than Ebh, such that the 
Universe has to be smaller than a critical size in order to hold the bound Equation ( 125 1 as is pointed in for the 
case of a vacuum dominated universe. 
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C. Type III Singularity 



This type of singularity is very similar to the Big Rip, in spite of the scale factor a{t) being finite at the singularity. 
The solution Equation (129) reproduces this singularity by taking — 1 < 7n < 0. The scale factor goes as: 



a{t) = Us exp 



Hi 



1 



{ts - t) 



where for simplicity we take Hq = 0. Then, for t — > tg, the scale factor a{t) 
singularity, let us write it in terms of the time instead of the scale factor: 



(144) 



— 7> fls. To see how Eq behaves near the 



(145) 



where m < 0. Hence, the Casimir energy diverges at the singularity, while Ebh oc a"~^ takes a finite value for the 
singularity time tg, so the bound is clearly violated long before the singularity. Then, in order to maintain the validity 
of the bound Equation (1251, one might assume, as in the Big Rip case, that GR is not valid near or at the bound. 



Even if quantum effects are included, as was pointed in J163' , for this type of singularity the energy density diverges 
more rapidly than in the classical case, so that the bound is also violated in the presence of quantum effects. 



D. Type IV Singularity 

For this singularity, the Hubble rate behaves as: 

H = Hi(t) + (ts - t)" H2{t) (146) 

Here Hi{t) and H2{t) are regular function and do not vanish at t = ts- The constant a is not integer and larger than 
1. Then the scale factor behaves as 

Ina(i)- / dtHi{t)+ / dt{ts-t)'^ H2{t) (147) 



Near t = ts, the first term dominates and all quantities like p, p, and a etc. are finite and therefore the bound 
Equation ( 125 1 would not be violated near the singularity. 



E. Big Bang Singularity 

When the matter with w > coupled with gravity and dominates, the scale factor behaves as 

a ~ t^^fA^ (148) 
Then there appears a singularity at i = 0, which may be a Big Bang singularity. Although the Big Bang singularity 



is not a future singularity, we may consider the bound Equation (125) when t ^ 0. Since n(l + w) > 2, the energy 
density behaves as p ^ a-"(i+™) and therefore the Casimir energy behaves as Ec ^ a^'"". On the other hand, we 
find Ebb ~ a"~^. Then when n > 2 or when n> 2 and w > 0, Ec dominates when a — 0, that is, when t 0, and 



the bound Equation ( 125 ) is violated. This tells us, as expected, that quantum effects become important in the early 
universe. 

Thus, in the above we have explored what happens near the future cosmic singularities. We have seen that in 
general, and with some very special exceptions in the case of Type H and Type IV, the bound will be violated if one 
assumes the validity of GR close to the singularity. Even if quantum corrections are assumed, it seems that the bound 
will be violated, although in the Big Rip case the singularity may be avoided when quantum effects are incorporated. 
It is natural to suggest, in accordance with Verlinde, that the bound on the Casimir energy means a finite range for 
the validity of the classical theory. When this kind of theory becomes saturated, some other new quantum gravity 



effects have to be taken into account. Hence, the universality of the bound Equation (1251 is not clear and may 
hold just for some specific cases, like the radiation dominated Universe. Furthermore, even in absence of a future 
singularity, a possibility allowed also for the case that u; < — 1 (see |192| 1. the bound may be violated, specially for 
the case of a "Little Rip" (see fl93j), a stage of the universe when the expansion would be strong enough to break 
some binding systems such as Solar System. That scenario, absence of singular points, can be very possible for a lot 
of DE models, and also in modified gravity (see |149l 11941 1195j ) 
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XI. CONCLUSIONS 

In this work we have considered both static sphericaUy symmetric solutions and massive charged and spinning 
solutions in f{R) theories of gravity in a vacuum situation within the metric formalism. The first case was analyzed 
in arbitrary dimensions, whereas for the second we focused on four dimensions. 

We first discussed the constant curvature case (including charged black holes in four dimensions) in static and 
spherically symmetric configurations. Then we studied the general case without imposing, a priori, the condition 
of constant curvature. For the uncharged case, we tackled this calculation through a perturbative analysis around 
the EH case, which made it possible to study those solutions being regular in the perturbative parameter. We have 
reproduced the explicit expressions up to second order for the metric coefficients, which only give rise to constant 
curvature (Schwarzschild AdS) solutions as in the EH case |1121 1113] . 

With regards to the massive, charged and spinning BH for f{R) gravity, we have derived the metric ten- 
sor that describes such objects. We presented how this solution differs from that found by Carter [119] 
by a proper redefinition of the spin, electric charge and vacuum scalar curvature and therefore that 
this type of solutions are supported in this kind of modified gravity theories. Further study of the 
Kerr-Newman metric allowed us to describe the different black-hole-type configurations derived from the presence 
of different horizons: usual, extremal and marginal extremal black holes, naked singularities and naked extremal 
singularities. 

On the other hand, we have also calculated thermodynamical quantities for the AdS as well as for KN configurations. 
We paid special attention to the issue of the stability of these kinds of solutions. As a very remarkable result, we 
have presented that the condition for a f{R) theory of gravity to support these kinds of black holes is given by 
^0 + f{Ro) < where Rq is the constant curvature of the AdS space-time. In fact, this condition also implies that 
the effective Newton constant is positive and consequently the graviton does not become a ghost. By employing the 
Euclidean action method, we revised how the mass, the energy or the entropy of these BH (AdS and KN) differ from 
those predicted in GR just by a multiplicative factor 1 -I- /'(Rq). This factor has to be positive in order to guarantee 
a positive mass and entropy for these kinds of BH. 

For the AdS configuration, we sketched the fact that the qualitative thermodynamic behavior of the AdS BH in 
f{R) gravity theories is the same as the one found by Hawking and Page. On the other hand, with regards to the KN 
configurations, the analysis of the BH heat capacity revealed how two different types of BH can be distinguished: fast 
and slow, showing the latter two phase transitions. With respect to the horizon structure in the KN configuration, 
we have shown how horizons can only exist for values of the spin lower than a maximum value a^ax^ and that from 
a certain positive value of the curvature onward, only above a minimum value amm- 

Then, for two paradigmatic f{R) models, we have investigated the stability of the different possible configurations 
that arise from the values of the free energy and the heat capacity, quantities that depend on the particular f{R) 
model under study. We have analyzed explicitly the rich thermodynamical phenomenology that characterizes these 
gravitational models with a complete set of different figures that were originally presented in |1121I113] and [1201 1121] . 
For doing so, we have studied the parameter regions in which BH in such models are locally stable and globally 
preferred. 

Furthermore, we have investigated the deep connection between FLRW equations and the first law of thermody- 
namics. The study of FLRW metrics is a very important issue as the universe seems to be well described by this kind 
of metrics, where eventually one should assume extra components to explain the cosmological evolution. However, as 
suggested above, modified gravity can even avoid the need to introduce dark energy, and even a kind of inflaton, to 
explain the whole cosmological evolution. Hence, FLRW metrics and their connection to thermodynamics turns out 
to be a fundamental issue. Then, as suggested in [5214M] . FLRW equations can be derived by assuming the entropy 
relation with the area of the apparent horizon. Even more, such a relation can be extended to /(i?) gravity when one 
assumes the modified relation for the entropy of the horizon. 

In addition, we have also reviewed the possible extension of generalized CV formula for multicomponent fluids, 
generalized in the sense that an inhomogeneous EoS (including modified f{R) gravity) was assumed. For some special 
cases the formula is reduced to the standard CV formula expressing the correspondence with 2d CFT theory. The 
dynamical entropy bound for all above cases was found. The universality of dynamical entropy bound near all four 
types of the future singularity, as well as the initial Big Bang singularity, was investigated. Note that the study 
of future singularities is a crucial step since a lot of dark energy models contain some kind of singularity. In fact, 
so-called viable modified gravities usually cross the phantom barrier {w < —1) (see jl96j ). and although it seems that 
future singularities do not occur, the "Little Rip" event may take place (see |151p . It was proved that except for some 
special cases of Type II and Type IV singularity, the dynamical entropy bound is violated near the singularity. Taking 
into account quantum effects of conformally invariant matter does not improve the situation. Hence, the dynamical 
entropy bound may not be universal and its violation simply indicates that the situation will be changed with the 
introduction of quantum gravity effects. 
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This deep connection between gravitation and thermodynamics in so many scenarios, as the ones reviewed here 
in extended gravitational theories and others (see |197j ). may suggest that gravity could be just an emergent force 
as suggested in [198] , which may imply the birth of new physics that may lead to the search of a more fundamental 
theory. 

Experimental checks to test the validity of a particular f{R) model can be applied in multiple realms of astrophysics. 
Among others, future data from CMB tensor perturbations |199j as well as gravitational collapse predictions |200l - 
I202| and the growth of cosmological perturbations (and subsequent structure formation) [17H5D] may constrain the 
available range of parameters for viable models ^203. .204] . In addition, the study of the weak lensing in f{R) gravity 
can also provide complementary constraints on the models and free parameters as pointed out in [205) . The general 
viability conditions for f{R) theories prior to model specifications are summarized in reference |114j . Furthermore, 
the study of star structure and evolution in modified gravity theories may lead additional and independent tests for 
this class of fourth order theories of gravity (see 206 ) . Finally, the aforementioned constraints may be enlarged not 
only by studying astrophysical BH stability but also if quantum gravity scale is accessible at particle accelerators. 
For example, it has been shown that LHC could produce about one microBH per second |207| I208j . Stability and 
thermodynamical properties of the produced BH might shed some light about the underlying theory of gravity. 
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